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Abstract 

We compute the cross section for leptoproduction of the real photon off the nucleon, which 
is sensitive to the deeply virtual Compton scattering amplitude with power accuracy. Our con¬ 
siderations go beyond the leading twist and involve the complete analysis in the twist-three ap¬ 
proximation. We discuss consequences of the target and lepton beam polarizations for accessing 
the generalized parton distributions from experimental measurements of the azimuthal angular 
dependence of the hnal state photon or nucleon. We introduce several sets of asymmetries, dehned 
as Fourier moments with respect to the azimuthal angle, which allow for a clear separation of the 
twist-two and -three sectors. Relying on a simple ansatz for the generalized parton distributions, 
we give quantitative estimates for azimuthal and spin asymmetries, discuss the uncertainties of 
these predictions brought in by radiative corrections, and compare them with experimental data 
as well as other theoretical expectations. Furthermore, we derive a general parametrization of the 
DVCS amplitudes in the region of small Bjorken variable. 
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1 Introduction 


The complicated long distance dynamics of hadron constituents, not controllable by well-dehned 
methods of perturbative QCD, was studied over decades by means of experimental measurements 
with inclusive lepton-hadron and hadron-hadron scattering. Information on the nucleon’s parton 
structure, originating from them, can be used to access inclusive properties of hadrons, like parton 
distributions. However, they are insufficient to constrain the detailed picture of the hadron wave 
function in full diversity of its manifestations. 

The parton distributions and correlation functions encode the long-distance dynamics of inclu¬ 
sive cross sections. The former accompany the interaction of a hard probe with a single parton, 
while the latter are involved in the rescattering off a multi-parton system. Both of them are given 
by the interference of light-cone wave functions with equal or different numbers of constituents 
in the amplitude and its conjugate, but the same total momentum flowing through each of them. 
This obviously simplihes the theoretical description on the one hand and, however, lacks a plethora 
of important features of strong interaction dynamics on the other one. 

As has been recently understood, a straightforward generalization of the aforementioned parton 
densities via new nonperturbative nucleon characteristics arises in exclusive two-photon processes 
in the generalized Bjorken region |^, e.g., in the Compton scattering with highly virtual incoming 
7 -quantum ^ ^ , and in the hard photoproduction of mesons . These are genuinely exclusive 
processes. Here one hnds off-forward matrix elements, as distinguished from the forward ones in 
inclusive reactions, of certain non-local quark and gluon composite operators, which give rise to the 
so-called generalized parton distribution (GPDs). They carry information on parton distributions, 
distribution amplitudes, and form factors and are of interest in their own right. Moreover, they 
can shed some light on the intriguing spin structure of the proton as they provide a way to get 
access to the parton’s angular momentum unreachable elsewhere P]. 

The deeply virtual Compton scattering (DVCS) is the hard photoproduction of a real photon, 
i.e., 7 *iV — '-yN'. Being a process involving a single hadron, it is one of the cleanest tools in the 
problem of constraining GPDs from data. The theoretical efforts and achievements are supported 
by pioneering experimental results from HERMES 0 , HERA ^ , and CLAS [§] collaborations, 
and encouraging future plans |^, |TT|, |T2 . 


This paper is a culmination of our previous studies initiated in Ref. |13| and completes the 
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theoretical background of the theory of DVCS by an exhaustive set of analytical results for the 
cross section to the twist-three accuracy for all possible hadron and lepton polarizations involved. 
Numerical studies performed here are done in the so-called Wandzura-Wilczek (WW) approxima¬ 
tion for the twist-three GPDs, developed by two of us in Ref. |^. Our present consideration is 
compatible with the one performed earlier in a fully numerical manner |^, see also Ref. |^. We 
also comment on the effect of multi-particle correlations, provided we cease using the WW rela¬ 
tions. The knowledge of power corrections to the leading twist-two results, with the twist-three 
terms being the hrst of them, is indispensable because the hard momentum of the incoming photon 
is not in the deep Euclidean domain for the present experiments. Therefore, the higher-twist con¬ 
tributions may affect the understanding of lowest order results. There is a number of theoretical 
uncertainties in phenomenological predictions for DVCS. Apart from the basically unconstrained 
shape of GPDs, which one models on the ground of previous theoretical considerations, these are 
higher-twist and -order effects, which alter the handbag approximation. The former include the 
dynamical and target mass corrections, see 0 for a partial result. While the latter have been 
computed in next-to-leading order (NLO) of perturbation theory in the strong coupling constant. 


see Refs. 


for coefficient functions and Refs. 


for evolution kernels. Their numer¬ 


ical signihcance was studied in Refs. ||2^, Q 0 . We will elaborate on these issues in the present 
study as well. 

Our presentation is organized as follows. The next section recapitulates basic results for the 
DVCS amplitude to the twist-three accuracy and introduces the WW relations for the twist-three 
quark GPDs. In section ^ we briefly discuss the gluon double helicity-flip distributions, vanishing 
in the forward scattering off spin-1/2 targets. Their peculiar properties propagate into a gen¬ 
uinely distinguishable angular dependence in physical observables. This allows to access them, 
uncontaminated by leading twist quark GPDs. The differential leptoproduction cross section is 
extensively analyzed with power accuracy in section ^ for all polarization options of the lepton 
beam and target. Based on these results, we discuss a way to unravel GPDs from cross sections 
and then turn to appropriate dehnitions of the most favorable physical observables, i.e., asym¬ 
metries. These asymmetries allow to extract separate components of the angular dependence of 
the cross section, and in this manner to project out distributions carrying information on the 
orbital momentum of constituents in the nucleon. In section p, we address the models for GPDs 
and discuss their diverse properties. We then study the magnitude of radiative corrections to the 
Gompton form factors (GFFs). Next in section |^, htting the model parameters to the recent HI 
data in the small Bjorken-variable region, we give phenomenological predictions for the HERMES 
and GLAS kinematics. Finally, we conclude. 
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2 DVCS amplitude in twist-three approximation 


Let us give a few basic definitions required for understanding of the Compton amplitude with 
power accuracy. Here we mainly recapitulate the results of our previous studies, where we have 
introduced the WW-approximation [|^ for the twist-three GPDs, and used it for preliminary 


studies of appropriate physical observables |]^ 


The DVCS hadronic tensor is given by the time-ordered product of the electromagnetic currents 
= e of quarks, having fractional charge Qi, which is sandwiched between hadronic 

states with different momenta. In leading order (LO) of perturbation theory it reads [H 


T,,{q,P,A) = -j dxP^-^{P2\Tj,{x/2)3,{-x/2)\P,) 
= —Pfiago-rPru^—— + (PP^Ppu + P) 


( 1 ) 


P V ] 

fip-*- a' <ju ) 


V, 


2p 


Pfia^^arqpP tu 


^1/ 


P ' U.U -u / uu ' I T~i ’ 

■ q P ■ q 

where we have kept all contributions up to twist-three accuracy. The three independent four- 
momenta are P = Pi -|- P 2 , A = P 2 — Pi, and q = {qi + q'2)/2, where the vectors Pi {qi) and P 2 (^ 2 ) 
refer to the incoming and outgoing proton (photon) momentum, respectively. The decomposition 
(|I|) is similar to the one used in deeply inelastic scattering. Indeed, the twist-two part of the 
generalized structure functions Vi and Ai corresponds to the conventional one Pi and gi (see the 
hrst paper of |^). The current conservation in the tensor decomposition (|^ is ensured by means 
of the projection operator 


Ppu g^u 


qipq2u 

qi ■ q2 


( 2 ) 


This is consistent with an explicit calculation of the amplitude (|l]) to twist-three accuracy ^ 


see also Refs. ||^, for spinless targets. 

In Eq. (|^ we have used conventional conditions on the kinematical invariants in the Bjorken 
limit, —q^ ^ P ■ q = large, = small, a generalized Bjorken variable ^ = —q^/q ■ P = hxed. If 
both photons are virtual, we would have an extra scaling variable rj = A ■ q/P ■ q, the skewe dn e.ss 
[|I|. The reality of the outgoing photon implies the presence of only one scaling variable namely, 
for q 2 = 0 one has 


h = -u f + 


A2 ' 


-1 


It is convenient to introduce variables closely related to the ones used in deeply inelastic scattering, 
namely, the incoming photon virtuality and conventional Bjorken variable 


Q^ = 




Xb = 




2Pi ■ qi 


( 3 ) 
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with the latter related to ^ via 


^ = xb 


1 + 


2Q2 


2- Xb + Xb^ 


( 4 ) 


In the above Eq. (|l|), V 2 p is expressed in terms of the vector Vip and axial-vector Aip fnnctions[|, 

( 5 ) 


v'2. = ■ n + 


2P-q^ 2P-q 

The amplitndes Vi and Ai depend on the scaling variable momentnm transfer A^, and hard 
momentnm of the probe Q^. The latter dependence is logarithmic and is governed by a generalized 
off-forward evolntion eqnation ||^, which is presently known to two-loop order [^, ^|. In 
order to simplify notations considerably, we will drop the dependence on and when it is 
not essential for the presentation and restore them in sections dealing with the modeling of GPDs 
and nnmerical estimates. 

A general decomposition of the vector and axial-vector amplitndes, in a complete basis of CFFs 
to twist-three accuracy, reads 


n = -Pp^w + Pp-^P + \ p 

q ■ P q ■ P ^ q - P 




p„ + + + (6) 


q ■ P 


q ■ P q ■ P 




q ■ P 
q ■ h 


P q ■ q-P^+^^P q-P 


q ■ P 

(7) 


where A^ = Ap -|- ^Pp, A^ = iep/^pq/P ■ q, and the Dirac bilinears are conventionally defined by 


hp = u{P2,S2hpU{Pi,Si), 
hp = U{P2,S2hp'y5U{PuSi), 


ep = UiP2,S2)tap^^U{P^,S^), 


A. 


2M 


U{P2,S2h^U{Pi,Si), 


( 8 ) 


with U being the nucleon bispinor normalized as U{Pi)U{Pi) = 2M. 

The CFFs introduced in Eqs. are given by a convolution of perturbatively calculable 

coefficient functions and a set of twist-two and -three GPDs via 

{h, 8, nl, si nls^} (0 = fdx x) {h, e, hI eI hIe^} {x, ^1=-^ , (9) 

{n,s,nl8lnls^_}{l = fdxC^+\^,x) {H,E,Hl,El,HlEl}{x,r]l=_^, 


where 1/Q^-power suppressed effects have been neglected, i.e., rj = —The GPDs {H ,..., E^} 
and {H,... ,E^} are related to the off-forward matrix elements of the vector and axial quark 
^ We use the conventions for Dirac and Lorentz tensors from Itzykson and Zuber e.g., _ _|_2 
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operators, see Eqs. (|116|) and (p.l7|) , respectively. We have implied that there is a summation 
on the right-hand side of the above equations over the quark species, so that they have to be 
understood as follows 


( 10 ) 


i: 

i=u,d,s 

with having perturbative expansion 

=do?+gdo + . 

To reduce the number of equations we have introduced unifying conventions F and for all 
twist-two and -three functions, which run over II, E, H, and E species as in Eq. At LO in 
the QCD coupling constant the coefficient functions read for the even (—) and odd (-I-) parity 


sectors 






doh -1 _ _ ,0 ^ 1 + x/J - *0 • 

According to |3, all twist-three GPDs are decomposed into the so-called WW term 
(later recalculated in Ref. |^) and a function F^^ that contains new dynamical information 
arising from antiquark-gluon-quark correlations: 


= F^^ + f: 




( 12 ) 


The WW parts are expressed solely in terms of the twist-two functions F = {H, E, H, E} and 
have the following form 


F 


WW 


= J 


-1 e 


W 4 


X y 


y 


d 


r^WW, 


dy 


X y 


dy 

d 


(y<0 (1_^2)(A2-Aij’ 


d 


f 0 = - -fW.\-,^-\\y^-i^\F{y,i)- 


1-1 e 




4M^F^{x,0 


(13) 


dy d^^ 

for the ‘-I-’ and ‘ —’ component of the twist-three GPDs in the WW-approximation, respectively. 
The minimum value of the momentum transfer squared is given below in Eq. (|31|). The 


functions F± specihcally appear for the spin-1/2 targets, cf. Refs. M, 1^, 
pion target, and read 


for the 


H^{x,^) = T 


A2 


X y 


X y 


4M2 J-i 


dy + E) {y, 0 - W:, ^ H{y, 0 


Ei{x,i) = ± j\yUw± 


Hi{x,i) = ±j^dyU[l- 


Ei{x,0 = F-^ j_^dy\W^ 


X y\ 

A2 
4M2 
'x y'' 


.re 

(H + E) {y,0-W^ 

W± 


X y 

X y\ ... A^ 


-y-\H{y,0 + 


- 

4M2 ^ 


e’e 

H{y,Oy 

'x y'' 






X y 


(H + E) (2/,0 . 


(14) 
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The hh-kernels introduced in these formulae are dehned by 


The antiquark-gluon-quark contributions, 


1 + y 


— 


" r l w f-l. ^SUy,u, -0 ± ^S-Ay,-u, -i) 


1-1 e J-i 


C U dy^ 


A' U 


can be read off from the parametrization of the corresponding operators and result in eight inde¬ 
pendent functions 


a± , C± . 






q ■ P ^ 


St±Aj^S±, 


q ■ P ^ 


(16) 


(see Ref. [jl^ for details). 

It turns out that only the difference T\—T^ enters the DVCS amplitude [Q, |^, |^, where, 

again similar to the previous dehnitions, the CFF P in Eq. stands for Ti, E, H, and S species. 
Therefore, only four new GPDs remain at twist-three level. In the WW-approximation, i.e., 
neglecting the antiquark-gluon-quark correlation, all the twist-three CFFs are entirely determined 
by the four twist-two GPDs if, E, if, and E. 


3 DVCS amplitude with gluon transversity 


Apart from the conventional vector and axial composite quark operators, whose matrix elements 
dehne the GPDs discussed in the previous section, there is a specihc gluon operator that requires 
the helicity of the hadron to be flipped by two units. It can show up in the scattering off a spin-1/2 
target due to a non-zero angular momentum of partons in the nucleon. The latter is a result of 
the non-vanishing i-channel momentum transfer A in the present settings. The corresponding 
operator forms the (1,1) representation of the Lorentz group and reads in terms of the gluon held 
strength tensors 

f^2) = G+pA2n)Til^.p„G„+{Kin ), (17) 


with the omitted gauge link. Here the totally symmetric and traceless tensor is 

1 


r, 




o (^ 9 up 9 i 


flp^PCT 


9p,(j9vp 9up9i 


1 p(J 


(18) 


7 






















The two-dimensional metric — n^n* — is defined in terms of two light-like vectors 

n and n *, such that = 0 and n ■ n* = 1. The off-forward matrix element of this operator 

is parametrized via four GPDs 


see 
dK 


3q, |39|| for an earlier discussion) 


A") = 4(P.n)-‘| 

?, A") ta-,e + f, A=) 


(19) 


+Et{x, A^)' 


P ■ q 

1 /7 ■ g 


2M2 


A^-rnp] - ET{x,ti,A^)^}U(P,), 


2M \P-q ^ ' '7 ' '2M 

The traceless symmetric projector r-*“ in Eq. (|T^ possesses the properties 


‘pa,p'a 


— '7~”^ — T~ 

' ocr:o'a' ’ uu’.oa ' i 


PxV]pG p(7]PxU ’ 


r-^ = n T-^ =2 

PxP]p(T ’ UV\LLU 


pV]pxU 


A simple calculation of one-loop diagrams |]38| , |3^ gives us the following result for the real final- 
state photon DVCS amplitude 


i=u,d,s ^ 
(+) 


»1 


( 20 ) 


into the 


with Tp = 1/2 and the coefficient function defined in Eq. (|^. Substituting Eq 
above expression gives the amplitude, which we will use in our computation of the cross section. 
We define the CFFs similarly to Eqs. @ via 

{nT,ST,nT,£T}iO = ^Tp /_'dxC'«(e,x){i7r,i^T,^T,^T}(x,0- (21) 

i=u,d,s ^ 

Below we will use unifying conventions for these CFFs, i.e., Tt = {^t, • • • 7 ^t}- 


4 Angular dependence of the cross section 


Now we are in a position to turn to physical observables, which give direct access to GPDs in a 
measurement of the five-fold cross section for the process e{k)h{Pi) —> e{k')h{P 2 )'y{q 2 ), 

da a^XBV 

dxBdyd\A‘^\d(j)dip 16 vr^ Q}\/\ 



This cross section depends on the Bjorken variable xb, the squared momentum transfer = 
(F 2 — Fi)^, the lepton energy fraction y = P\- qi/Pi ■ k, with qi = k — k', and, in general, two 
azimuthal angles. We use throughout our presentation the convention 


e = 2x 


M 
'"‘Q ■ 
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Figure 1: The kinematics of the leptoproduction in the target rest frame. The z-direction is 
chosen counter-along the three-momentum of the incoming virtual photon. The lepton three- 
momenta form the lepton scattering plane, while the recoiled proton and outgoing real photon 
dehne the hadron scattering plane. In this reference system the azimuthal angle of the scattered 
lepton is (pi = 0, while the azimuthal angle between the lepton plane and the recoiled proton 
momentum is (pN = (p. When the hadron is transversely polarized (in this reference frame) 
S± = (0, cos^, sin^, 0), the angle between the polarization vector and the scattered hadron is 
denoted as (p = — (pjsi. 

In Eq. d^), (p = (pN — <pi is the angle between the lepton and hadron scattering planes and 
ip = $ — (pN is the difference of the azimuthal angle ‘P of the transverse part of the nucleon 
polarization vector S, i.e., S_\_ = (0, cos^, sin^, 0), and the azimuthal angle 07v of the recoiled 
hadron. Our frame is rotated with respect to the laboratory one in such a way that the virtual 
photon four-momentum has no transverse components, see Fig. |l[ We £x our kinematics by 
choosing the z-component of the virtual photon momentum to be negative and the positive x- 
component of the incoming lepton: k = (E, E sin 6*;, 0, E cos 6*;), qi = (o'?, 0, 0, —|gf |). Other 
vectors are Pi = (M, 0,0,0) and P 2 = (E 2 , IP 2 I cos^sin^Ar, IP 2 I sin^sin^Tv, |T* 2 | cos6*Ar). The 
longitudinal part of the polarization vector is S'lp = (0, 0, 0, A). 

The amplitude T is the sum of the DVCS Tbvcs and Bethe-Heitler (BH) Tbh amplitudes. The 
latter one is real (to the lowest order in the QED hne structure constant) and is parametrized in 
terms of electromagnetic form factors, which we assume to be known from other measurements. 
The azimuthal angular dependence of each of the three terms in 

= |%h|^ + l^vcsP + (23) 
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with the interference term 


^ — ^VCS^*H + ^ 


-dvcs^bh , 


Tr 


(24) 


arises from the contraction of leptonic and handronic tensors (see also [^). In onr frame these 
contractions yield finite snms of Fonrier harmonics, whose maximal freqnencies are defined by the 
the rank-m of the corresponding leptonic tensor in the incoming lepton momentnm k^. Note, 
however, that in the polarized part of the leptonic tensors, proportional to A times the e-tensor, 
one has one fonr-vector less than in the nnpolarized part. Thns, the highest harmonic which is 
proportional to A will be cos / sin([m — !]<;/)) instead of cos / sin(m0). The parity and time reversal 
invariance provide further constraints on the Fourier coefficients. 

The BH term ITbhP, squared DVCS amplitude l^bvcsPi and interference term X read 


\T |2 — 


a;|2/2(i + e2)2A2X’i(0)'P2(0) 


Co “ + £ cos {ncf)) + sin (0) 


n=l 


\Tt 


DVCS — 


1 = 


y2Q2 

±e® 


c, 


DVCS 

0 


+ !:[■ 


n=l 


DVCS 


cos(n0) -F sin(n0) 


Xb2/^A2X’i(0)X2(0) 


+ [cncos(n0) + s^sin(n0) 


n=l 


(25) 

(26) 
(27) 


where the -|- (—) sign in the interference stands for the negatively (positively) charged lepton 
beamQ. The results for the Fourier coefficients, presented below, show that the generation of 
new harmonics in the azimuthal angular dependence is terminated at the twist-three level. The 
coefficients cf, s\ as well as arise at the twist-two level, and their dependence on GPDs has 

been elaborated in Refs. ^|. The rest provides an additional angular dependence and is given 

, cf, and sf, GPDs. The harmonics 


m 


terms of twist-two, i.e., cj, and twist-three, i.e., c 


DVCS „DVCS 


1 


’1 


proportional to cos (30) [cos (20)] or sin (30) [sin (20)] in the interference [squared DVGS] term 
stem from the twist-two double helicity-flip gluonic GPDs alone. They are not contaminated 
by twist-two quark amplitudes, however, will be affected by twist-four power corrections [^ . 
We neglect in our consequent considerations the effects of dynamical higher-twist (larger than 
three) contributions. They will give power-suppressed corrections to the Fourier coefficients, we 
discussed. 

There is an important difference between the interference term and the squared DVGS ampli¬ 
tude. The former has a contaminating 0-dependence due to the lepton BH propagators. 


= {k- q^Y = Q^ + 2k-A, Q^V2 = {k- A)^ = -2fc ■ A + A^ 


(28) 


^Note that only in the case of massless leptons, a partial cancellation of propagators occurs in the squared BH 
term |^ . 
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where 


k-A = 






(29) 


The 1/ Q-power suppressed kinematical factor K appearing here also shows up in the Fourier series 


4xb(1 - xb) + 


-b) (i - y - A) (i _ 


S" 


(30) 

(with the plus sign taken for the square root in Eq. (^9D ). It vanishes at the kinematical boundary 
A^ = ATj^, determined by the minimal value 


_ ^ P(1 - ^b) (1 - vrr?) + 

““ 4 xB(l-XB) + e 2 


M^xl 


1 - xb + xbM ^/Q 


2 ’ 


as well as at 


y{x, Q^) = ?/max = 2 


v^TT^ - 1 ^ M^xl 

-^ ^ I — 


e- Q2 

The square of the transverse momentum transfer is given by 


(31) 


Al«(l-5^)(A^-Ai„J, 


up to corrections suppressed by the hard-photon virtuality. 

According to Eqs. (^) and we introduce the following parametrization: 


Vi = 


y{l + e^) 


{J -|- 2K cos(0)} , 


A2 

^2 = 1 + — + 


Q2 1/(1 + 62) 


{J + 2K cos(0)} , (32) 


where 


As we see, the denominator of the u-channel lepton propagator, i.e., Vi, can be of order 1/Q^ at 
large y. In the Bjorken limit it behaves like {1 — y). Moreover, if the outgoing photon is collinear 
to the incoming lepton, it vanishes. Of course, the photon then lies in the lepton scattering plane, 
i.e., = 0 + TT = 0, and both polar angles coincide with each other. The latter condition is 

fulhlled if 


_ Q2 + a2 
y - ycoi - Q2^ ^^2 


A2 

1 + (1-xb)-^, 
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(with j/coi < I/max) and this leads to the equality 


'^\y=yco\ 


2K\ 

^-f^|y=ycol 


( 33 ) 


2(l-a:B)A2 



Vl + + 


4xb(1 - xb) + 


4(1 - xb) 


Q2 


Obviously, for large y the squared BH (^) and interference (^) terms are enhanced with respect 
to the squared DVCS one (PB|). Furthermore, the expansion of Vi in Q is not justihed, and, thus, 
the Fourier analysis of experimental data must be modihed. For small y it is legitimate to expand 
Vi and 1^2 in power series with respect to 1/Q. This generates higher harmonics suppressed by 
powers of K. 

The computation of the Fourier coefficients in Eq. (p5|-p7|) is tedious, however, relatively 
straightforward (see Refs. ^ ^]). After an extensive algebra we come to the Fourier har¬ 
monics summarized in the following sections. We will employ the conventional definition of the 
lepton helicity, i.e., A = -|-1 if the spin is aligned with the direction of the lepton three-momentum. 
To have a compact notation, we write the cross section for a polarized target as 


da = diJunp + cos(6') dcrLp(A) -|- sin(6*) daTB{ip ), (34) 

where the polar angle 9 appears in the decomposition of the spin vector S = cos(6')S'lp(A) -|- 
sin(6')S'x($). The same decomposition will be used for the Fourier coefficients. 


4.1 Bethe-Heitler amplitude squared 

This part of the leptoproduction cross section is expressed solely in terms of Fi(A^) and F 2 (A^), 
the known Dirac and Pauli form factors of the nucleon. For diverse target polarizations the Fourier 
coefficients are: 

• Unpolarized target: 


..BH 


0,™p - 8A'^ I (2 + 3P) A + 24 (F, 

+ {2-yf\ {2 + e^) 


+ Ao 


AxlM^ f A^\ W A2 

1 + ) + 4(1 — Xb) ( 1 + Xb^^ 


(35) 


A2 


+ Axi 


a:B + 1 - Xb + 


1 - 


A 

Q 


Q\ V 

2\2 ^4 

-Xb(1 -2xb)7^ 




2- , . 
(i^i + F2f 


. ^ 4M2^2 


+ 8(l + e2) (l-y- 


eV' 


2eM 1 - 


A2 ' 
4M2 


^2 - ^F^' 
. ^ 4M2^A 
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cf,S„p = 8A'(2-,)|(±|^-2xB-^'|l^’^ 


2 1 / I?2 p2 

4M2 2 ^ 


(36) 


A2' 


+ 2xb ( 1 - (1 - 2xb)-^ j {Fi + F2f 


c?,“„p = 8 x|A- (a^ - ^A|) + 2 (Ap + Ap)^ 


(37) 


Longitudinally polarized target: 


cjj!p = 8 AAp:b(2 - y)v^^F^ (A, + App) (1 

4M2 I' ^ 


2 I Q2) 4jvf2 


2 - xb 


A2 / A2\ A^' 

- 2(1 - xb)^^ + 1 ^ - 2^")^ 


{Fi + F 2 ) 




xlM^ /. A2\" ,, , /, A2' 


A2 ( 1 + xb-^ 



c™p = -SXAxbvK 


A2 

4M2 


(i^l + i"2) 


A2 


A2' 


2M2 r Q2 


A^' 


(38) 


1 — Xb + 1 (-^1 + -^ 2 ) 


1 + xb - (3 - 2xb) ^1 + xb^^ - ^2 


Transversely polarized target: 


A2\ 4xlM^ 


A2 ^ 


(39) 


„bh _ 

^0,TP ~ 


Q Al + e'^K 

-8A cos {cp){2 - y)y4- / ^ (F^ + F^) 




XrM2 A2\ _ Xb 




A2 I ^ Q2 ) + 2 


„bh _ 

'-l,TP ~ 


2iL2Q2 


/ e2'y2 M I 

-I 6 ACOS (^)ib!/\/i - !/ - —-g Vr+^(A, + -Ri) 1^2 j-j _ ^ 

j^aJ + (i + t=)iB (i - |j) (a. + I. 


(40) 


XB (1 - |j) Ap 


gSH = i6y sin {ip)yxl\l 1-y- 


e2|/2 M 

~T~q 


+ <t(i-|))(Ap + A.)(ap + ^Ap' 


(41) 

(42) 


4.2 DVCS amplitude squared 

|7bvcsP is bilinear in the CFFs, and its coefficients read in terms of functions, which are 


specihed in section WA 
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Unpolarized target: 


^DVCS 

"'OjUnp 


= 2{2-2y + y^)C^l^^ 


^DVCS 'i 
-'ljUnp I _ 

^DVCS I ^ 


%K 


’ l,unp 
^DVCS 

'^2,unp 


2-Xb 

4Q^K^ 


2 - 1 /' 

-Xy 


3fJe 


9m 


C, 


DVCS 

unp 




M2(2 -xb) 
Longitudinally polarized target: 




(43) 

(44) 

(45) 


^DVCS 

'^0,LP 


'„DVCS■ 
'-1,LP 

„DVCS 
d’l.LP , 


„DVCS 

'’2,LP 


2XAy{2-y)CE^^^ 



^DVCS 

'-'LP 


4AQ^iL^ 9m (T T* 

^muj.Lp 

x-b) 




Af2(2 


(46) 

(47) 

(48) 


• Transversely polarized target: 


„DVCS 

-0,TP 


' mvcs 

A,TP 

„DVCS 
. ■’1,TP 


QK 


-Xyi2 - y) cos((y9)CT^^® + {2-2y + y^) sin((p)9mCT^^“ J^*) , 


DVCS 


' „DVCS 
C2 TP 

„DVCS 
. *2,TP 


My/T^ L 
AQK"^ 


X 


M(2 - xb)a/i - y 

-Xy ] f 3 
2 - 2 / 


(49) 

(50) 


COs(:p) 


Ctp+® + sin((p) 


9m 


2 - 2 / 

Xy . 


9m 

3fJe 


C 


DVCS 

TP- 




AQy/T^K 

M{2 - xb) 


9m 


sin((p) 

cos((p) C°Tp+ 




(51) 


The harmonic is given in terms of twist-two CFFs JF = {H,S,H,£}, defined in Eq. (^. 

The twist-three coefficients and arise from the interference of twist-two CFFs with 

‘effective’ twist-three ones, 


jreff ^ _2^ 


1 + e 




(52) 


where 2F|. are defined in Eqs. (|^-]T). These Fourier harmonics have the same functional depen¬ 
dence on CFFs as the leading twist-two ones [Q. However, this is not the case for the Fourier 
coefficients and induced by the gluon transversity CFFs (pi]) . 
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4.3 Interference of Bethe-Heitler and DVCS amplitudes 

For the phenomenology of GPDs, X is the most interesting quantity since it is linear in CFFs. 
This simplihes their disentanglement from experimental measurements. The Fourier harmonics 
have the form: 

• Unpolarized target: 


( 2 -!/)" 


4 ,„„p = - 8(2 - (F) + |^(1 -y)(2- xb) (cJ., + ACj„„) (T) \ , 


I::} ^ 


_ iwr" (-(2-9)1 ^ 


^2,unp 


53m 


2-xb[ \y 

r 8Q2iF3 

Up = 


c„\p {r”') , 


(53) 

(54) 

(55) 

(56) 


Longitudinally polarized target: 


Cg^LP = —SXAy^e 


{( 


1 -y 


A 


+ 2 ) K^CIb (^) + -^(1 -y){2- xb) {cIb + ACIb) (^) \, 


'n,Lp'l i ~ 

aJ *“12-29 + 9^ 


3fJe ' 
Q'm 


cIp (j^) , 


c? 


r = 

’2,LP ) 


- 

Oq T U - 


16AX:^ r -\y 
2 -xb [2-y 
SAQ^K^ 


I cl, {F”) . 

Am 8 / 


3-LP M^{2-XBy 
Transversely polarized target: 


^mCj^Lp {J't) • 


(57) 

(58) 

(59) 

(60) 


Cq.tp = ~ 


Q 


- \y cos((^) 3fJe| + 2^ C|p+ {T) + AC|p+ (J^) 

+ {2 — y) sin((/i) Am| — —^C|p_ (JF) -h AC|p_ (JF) 
11-1/ 


4, 


TP 


’1,TP , 


8M^/l - y 

Q 

, , f -\y{2 - //) 1 f 3fJe 1 ^ r o _ o., , .,2 

X Icosfo?)-! n r ■) r 1^TP4- + sinfw) , ^ , , 

^^^\2-2|/ + |/2/\AmJ ^^>\\y{2-y) 


2 — 2y + y^\ j Am 


(61) 


'-'TP- 


(62) 

A) , 
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cf TP 1 16M^1 — yK 

/ Q(2 Xb) 


X 


003 ( 99 ) 


-\y' 

2 - 1 /, 

8Qv/r^ir2 


3?e 




Ctp+ + sin(99) 


2 - 1 / 

A// 


'^S.TP — ^^^2 COs(99)^5mCT^TP+ (-^r) , 


M(2 

8Qv/r^ir2 

3,TP ^(^2 - Xb)^ 


Co rpp - 


sin(99)9mCTTP- (^ t ) 


Q'm 



( 63 ) 


(64) 

(65) 


The higher twist-three harmonics, i.e., cf and sf have again the same functional dependence 
as the twist-two ones. However, this is not the case for cj, which depends only on twist-two CFFs 
JF, and for cf and sf, induced by jFp CFFs. 


4.4 Angular harmonics in terms of GPDs 

The Fourier coefficients displayed above are expressed in terms of the coefficients C. They depend 
on GPDs, integrated over the momentum fraction, and are functions of the kinematical variables 
xb, A^, and Q^. For the harmonics involving H, E, H, and F'-type GPDs they are: 

• Squared DVCS amplitude: 

^*) = ~ + ™*) - 4+ sn*^ 


Clp (F, E*) 


1 

{2-XBy 


4(1 - Xb) (nn* + nn*) - x| (ns* + sn* + ns* + sn*^ 

-XB(f + (2 -xb)^) (££• + ££•)}, 


(67) 


^DVCS /-77 79*'l 
C"pp+ J 


^DVCS 

Utp- 


{E,E*) 


1 

(2 - xb)2 
2 

(2 - xbY 


2 xB(ns* + sn*) - 2{2 - XB)(ns* + n*s) + xliss* + ss*)j, 

(2 - XB){ns* - sn*) - XB{ns* - sn*)]. (68) 


• Interference of Bethe-Heitler and DVCS amplitudes: 
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Here a part of the result at the twist-two level is expressed in terms of the functions, which show 
up in the lowest twist approximation, and they read 




- 

'-'LP — 


+ F,) (« + f f) + (f F, + g , 


Ctp-i- ~ (-^1 + -^ 2 ) 


^2 

Xp 




(n + 

2 - XB \ 2 


+ ^(4^^F2?{- fxBF,+ ^ 


- 

'-'TP- — 


4M2 1 2 - XB 
1 


2 - xb' 


-Fo\£ 


A2 


2 - xb 


XpFi - (1 -XB)77yF2 \ n + 


M2' 


r A2 

\4M2 


Xt 


-FA£ 


^2 

Xp 


-(F 1 + F 2 ) + 


A^ 


(2 — Xb)Fi -|- 


-£ 


xt 


2 - xp' 


2-xp^"J" 2-xp'^^ ' V’” ' 4M2 ^ 

The addenda arising in power-suppressed contributions are dehned as 


AC^ = ■ 

""P 2 - Xp 


{F^ + F2)\-^{n + £) + n\, 


2 - Xp 


AC^P = 


AC|p+ = 


+ F 2 ) {h + ^£ + 

2 — Xb I 2 2 — Xb 


H + ^g 


'f,W- 


M2 


, Fi + —F 2 ]£} , 

2-xb\ 2 ' ' ^ ’ 


A2 


^^TP- = {F 2 H - F,£) . 

Let us now list the coefficients involving the gluon transversity: 
• Squared D’VCS amplitude: 


(69) 

(70) 


(71) 


(72) 

(73) 

(74) 

(75) 


^DVCS 

'-'T,unp 


(2-xb)2 


Ht {2 - xb)£* - xb£* -2{2-xb)'Ht 


A2 


^DVCS 

'-'T,LP 


(2-xb)2 


—£t 

|?^T 

—£t 

F£t 


(2 - xb)'H* - xb'H*\ + £t xb{FC + £*) - (2 - xb)^^ 

(2 — xb)T* — xbT* -|- 'hij' 


(76) 


2(2 - XB)7f* - Xb I^xb - £* 


(77) 


XB^f* - (2 - XB)7f* + -^{£*+ £*) 


(2 - xb) ( n* + - XB ( 77 * + 


^DVCS 

'-'T,TP+ 


(2 - xb)2 


I 4(1 — xb)FCt ~ x\£t + xb( 2 — XB)7^r {fL* + Fi*^ 


(78) 
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-Ht {x"^ + (1 - 


-Xb 

-\-£t 


a;B^T + ( y + (2-xb)^ 


Si 


{£* + £* 


(2-xb) + 


2 Yb 

B o 


A2 ' 
4M2 




^DVCS 
TP_ 


(2-a;B)2 


I 4(1 — xb)'Ht ~ x^St + xb{‘2 — xb)£t 


-2nT (^xI + {1-xb)^'^ {HYsn 
-Xb (^xbHt - (2 - (S* + £*) 

Xg + (2 - Xb) 


-£ 


A2 \ A2 ~ 

^ ^£* + xl-^£* 


4M2 


4M2 


Interference of Bethe-Heitler and DVCS amplitudes: 


( 79 ) 


C; 


T,unp 
■I 


— —F2'Ht + 2 ( + 


Xb 


A2 

4ii^ 


F2 ) SCt + FiSx , 


^T,LP — -^2 ^ Ti-T + + £t + £t) r + Fi^Xb'Ht + £t) 1 


X 


Ct,tp+ — (2-^1 + xbF2)'Ht + a^B ( — — 


A2 


2 4M2 


xt 


+ < xb-Fi + — —1“ (2 — Xb) 


A2 \ 
4M2 j 
A2 


2(-Fi + F2)'Ht + -^2^-7 


-^2 ^ , 


Ct,TP- ~ (2-^1 + XbF2)'Ht — (2 — ^b) ^^2 + F2£t + (^B ( -^1 + ) £t ■ 


A' 


(80) 

(81) 

(82) 


'4M2 


(83) 


This set of formulae is the complete result for the real-photon leptoproduction cross section in 
the twist-three approximation. Below, presenting quantitative estimates, we will not discuss the 
case of transversely polarized target, therefore, the integration with respect to (p gives 27i on the 
right-hand side of Eq. 


4.5 Structure of effective twist-three GPDs 


Note that the twist-three GPDs, having generic discontinuities at |x| = ^ [^, |^, enter the CFFs 
in a singularity-free combination (^). Using the integrals of the type 


"i dx 


-lU 


X y 


sign(r7) 


In 


V + ^ 


-i\ri\^ — x — i0 \V^ V J VFy — V — 
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in Eq. the convolution of the coefficient functions and the fh-kernels results into an 

explicit form for the effective twist-three function, 

= —(^(0 + 2^3^ dx C-3WK. x)F(x. 0 H- (84) 

-2( j^^du J^^dxC‘'°‘’((,x,u) (SJ)-!, -u, -J) - Sp{x,u, -{)) |, 

where T is dehned in Eq. (|^). The summation over quark flavors is implied, as in Eq. (pT]), and 
we have used new notations for the coefficient functions weighted by the quark charge 


C 


3(t) 

(0)i 





Q\ 


In ■ 


2^ 


^ + X ^ — X — iO 
d'^ I + u 


^ + X — iO' 


Q 


dx'^ ^ + X 


In 



(85) 

( 86 ) 


The — (+) sign stands here for the vector (axial) sectors, i.e., for {7-f(7-f), S{S)}. Finally, the 
functions JF-*- are 


W^(0= f_^dx{iC^l-Hx,i)(H + E)(x,i)-C^^*'>(x,i)H(x,i)] , (87) 

£^(0 = f_^dx{(,C’^'--Hx,C)(H + E){x,i)-C^^*\x,QH(x,(,)} , 


^^(0 = / dxU 1 - 




A2 ' 
4M2 


(x, i)H{x, 0 + 0 {H + E) (x, 0 


^^(0 = 


dx {C3(+)(x, Off (x, 0 - e (a:, 0 (H + E) (x, 0} 


We should note that the kinematical factor (A^ — A^;^)“^ in Eq. (|8^) cancels out in the hnal 
results for the angular harmonics calculated in sections |4.2| and [4.3| . Therefore, it does not lead 
to a singular behavior of the Fourier coefficients in the cross section (^). For all interference as 
well as for unpolarized and longitudinally polarized squared DVCS terms this cancellation already 
appears in the C functions given in sections see Ref. |1^ for explicit examples. However, in 
c?TP^ and fhis kinematical singularity is annihilated by the prefactor in Eq. 


5 Physical observables and access to GPDs 

As we have seen, the cross section for leptoproduction of the real photon possesses a quite rich 
angular structure. The goal of experimental measurements is to pin down the GPDs and this 
requires a clean disentanglement of different components of the cross section (p2D. A general 
strategy, which follows from the analytical results, given above, will be discussed in section o 
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sector 


harmonics in X 

extraction 

P of 

A(j_ behavior 

twist 

C’s 

unp 

LP 

TP. 

TP, 

of CFFs 

Q-^ 

unp, LP 

TP 

two 

3f?eC(J^), AC(J^) 

Cl ; CQ 

Cl 5 Cq 

Cl, Co 

Si, - 

over compl. 

1,2 

1,0 

0,1 


AmC(J^), AC(J^) 

Si, - 

Si, - 

Si, - 

Cl, Co 

over compl. 

1,2 

1,0 

0,1 

three 

3fJeC(J^"®) 

C2 

C2 

C2 

S2 

complete 

2 

2 

1 


AmC(J^"®) 

S2 

S2 

S2 

C2 

complete 

2 

2 

1 

two 

3f?eCT(XT) 

C3 

- 

- 

- 

1 X 3fJe of 4 

1 

3 

2 


AmCr(XT) 

- 

S3 

S3 

C3 

3 X Am of 4 

1 

3 

2 


Table 1: Fourier coefficients cf and sf of the interference term defined in section |4.3|, while the 


corresponding C coefficients are given in section 4.4 


separation of twist-two and -three sectors^], measuring the CFFs, and access to the GPDs. To 


go along this line, we introduce appropriate asymmetries in section |5.2| and demonstrate how a 
Fourier transform can distinguish between the interference and squared DVCS contributions. 


5.1 Unraveling GPDs 

So far we have introduced eight CFFs at the twist-two level, with four of them from the gluonic 
transversity contribution. Another four new CFFs arise at the twist-three level. These sectors 
can be separated due to their characteristic azimuthal dependence as summarized in Table ^ for 
the interference X and in Table ^ for the squared DVCS amplitude, respectively. 

Let us first examine the issue of the dominance of each of the three terms (p5- 0) in the 
leptoproduction cross section (P^) in different kinematical regions. To do so we have to know 


the functional dependence of the Fourier coefficients given in Eqs. (p5| - p7D on scaling variables 
and transferred momenta. For instance, apart from an explicit XB-dependence of multiplicative 
prefactors there is one hidden in CFFs. It will be argued in section ^ that the latter behave like 
in the small-XB region. Thus, for general kinematical settings we expect for the scattering on 

XBcf/iF. Taking now 


, and (H^) that q 


..BH 


~ XrC, 


2 mvcs 


B'-O 


the unpolarized target from Eqs. 
the kinematical prefactors in Eqs. (p5| - p7D into account, and the behavior of the BH-propagators 
(|32D , we realize that the ratio of the DVCS to BH amplitude behaves like (1 — y) jQ}. 
Obviously, for small (large) y the DVCS (BH) one dominates. As compared to the squared 

^We talk about the twist-two and twist-three sectors, however, beyond the twist-three approximation each of 
the associated observables gets corrected by higher-twist, i.e., power-suppressed contributions. 
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Co 
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Cl 
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Si 

4 X 3?e/3m 
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Si 

Si 
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Cl 
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1 
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3f?eCT(J^T, 
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- 

1 X 3?e/3m 

2 

2 

- 



- 

S2 

S2 

C2 

3 X 3m/3?e 

2 

2 
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Table 2: Fourier coefficients and of the squared DVCS amplitude defined 


in section 4.2, while the corresponding C coefficients are given in section 


amplitudes, the interference term has an additional factor y'A^/A^. Note that the beam spin- 
flip contributions provide always an additional damping by the factor y. For the unpolarized 
or longitudinally polarized target, higher harmonics in any of the three terms are suppressed 
by powers of K. However, in the case of gluonic transversity, this goes hand in hand with an 
enhancement by It is important that lower harmonics in the interference term, i.e., cj, 

appear at the twist-three level. Since, they are not proportional to the factor K as compared 
to Cl, which is, they can be rather important and sizable close to the kinematical boundaries. 
Fortunately, they only depend on twist-two CFFs, thus, we listed them in Table ^ in the twist- 
two sector. In the case of transversely polarized target, we observe that both higher and lower 
twist-three harmonics are suppressed by one power of K in the interference term. 

The analytical structure and simple counting rules, we gave, provide a guideline to separate the 
three different parts in the leptoproduction cross section. In single spin-flip experiments, which 
give access to the imaginary part of CFFs, the BH cross section drops out, while in unpolarized 
or double spin-flip experiments it does not and has to be subtracted. This can certainly be done 
for not too small values of y. The interference and squared DVCS terms have different azimuthal 
angular dependencies due to the presence of the BH-propagators in the former. In principle, this 
fact can be used to separate them by a Fourier analysis. However, this method requires from the 
experimental side very high-statistics data and from the theoretical side a better understanding 
of twist-four contributions. 

Due to the different charge parity of individual components, it is possible to use the charge 
asymmetry to separate the interference and squared DVCS terms. The interference term is charge- 
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odd and can be extracted in facilities that possess both positively and negatively charged leptons 


4t3|, i.e., 


da~^ — da oc 


XBy^VMV2{(j))A^ 


cf + 


C5 

[cn cos(n0) + si Sin(n0) 


n=l 


( 88 ) 


Its measurement and consequent extraction of separate harmonics provides the real (unpolarized 
or double spin-flip experiments) and imaginary (single spin-flip experiments) part of linear com¬ 
binations of twist-two and -three CFFs. An explicit projection procedure of these harmonics will 


be discussed below in section |5.2.1| . Moreover, the charge-even part is given by the sum of the BH 
and DVCS cross section. The subtraction of the former gives then the Fourier coefficients of the 
latter: 

1 ' " 


d'^a + d a — 2d^^a oc 


..DVCS 

''0 


+ 


El' 


.DVCS 


cos(?7,0) -1- sin(?7,0) 


(89) 


In the case of the four CFFs T = {Ti, T, Ti, T} we have eight observables given by the first har¬ 
monics cos(0) and sin(0) of the interference term, which are accessible away from the kinematical 
boundaries in polarized beam and target experiments. Thus, experiments with both longitudinally 
and transversely polarized target can measure all eight Fourier coefficients cf yj and si and, thus, 
also 3?e/AmC^ with A = {unp, LP, TP^,, TP^,}. Knowing these C functions, we can invert them to 
obtain the CFFs: 
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Consequently, the four Fourier coefficients as well as the four twist-two DVCS coefficients 


can serve as experimental consistency checks. Alternatively, they can be used to extract 
CFFs. Thus, experiments with longitudinally polarized target have the potential to extract the 
real part of all four CFFs as well as two linear combinations of their imaginary parts from the 
interference term alone. The missing two imaginary parts could then, in principle, be obtained 
from the DVCS cross section, i.e., by measuring and 

For a polarized beam and target with all polarization options, the real and imaginary part 
of all four CFFs in the twist-three sector can be extracted from the interference term alone by 
projection onto the cos(20) and sin(20) harmonics and using Eqs. ( p^ - p3| ). Alternatively, knowing 
the twist-two sector and having only a longitudinally polarized target, one can employ in addition 
the squared DVCS term, i.e., its cos(0) and sin(0) harmonics, to access the full twist-three sector. 

For gluonic transversity, the cos(3(^) and sin(30) harmonics in the interference term can only 
provide one imaginary and three real parts of certain linear combinations of Missing informa¬ 
tion can, in principle, be obtained from the cos(20) and sin(2(;/)) harmonics of the squared DVCS 
term. Note that here again a polarized beam and target with all polarizations is necessary. More¬ 
over, the gluonic transversity is suppressed by Os/tt, so one expects a stronger contamination by 
twist-four effects H- 


As we discussed, a combination of the charge asymmetry with different nucleon/lepton polar¬ 
izations and projection of the corresponding harmonics provides, at least in principle, a way to 
explore the real and imaginary part of all CFFs (cf. P, ^). This gives maximal access to all 
GPDs, which enter in a convolution with the real or imaginary part of the coefficient functions, as 
we have established above in Eqs. and (pT]). Since these formulae cannot be deconvoluted in 
practice |^, one has to rely on models with a set of free parameters, which has to be adjusted to 
experimental data on CFFs. This issue will be discussed below for different available experiments. 


5.2 Asymmetries 

The measurements of the cross section (^) in different setups, as discussed in the preceding section, 
would directly lead to determination of the CFFs. However, on the experimental side it is simpler 
to measure asymmetries since they escape the normalization issue. Thus, we now discuss the 
separation of twist-two and -three sectors in their terms. The charge and lepton-spin asymmetries 


used in our previous studies, see, e.g., |^, were restricted to the leading twist approximation for 
the amplitudes and cannot serve the purpose since they are inevitably contaminated by power- 
suppressed effects. Due to our poor knowledge of multi-particle correlations inside the nucleon. 
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this affects the theoretical predictions in an nncontrollable manner 
The charge asymmetry 




/■3V2 _ ^-^unp' 

U —#— 


L ^ - d4 —■ 


for nnpolarized settings contains the contribution of all harmonics due to the presence of the 
non-negligible 0-dependence of BH propagators: 


with 


y„(2A7J) oc 


« i: y„(2A7 j)< 

n=0 


/■^/2 COs(?T.0) /•37r/2 

l-w/2 Pi ( 0)^2 ( 0 ) A /2 


d4>- 


cos(n0) 


(96) 


(96) 


'PiW'P2{<t>), 

while the normalization is not affected by twist-three corrections. If the final photon is collinear 
to the incoming (massless) lepton, Pi(0) is peaked at 0 = vr. Thus, the ratio Iin/Iii approaches 
plus or minus 1 and all harmonics contribute on equal footing. However, since K is then of order 
Q^, only cj and give essential contributions at the same order in A^/ Q^. Thus, for this 
asymmetry cf may gives an essential effect for large y, since the twist-two part becomes small. 
For y -C ycoi, all twist-three harmonics are suppressed in addition by the iF-factor. From the 
expansion of „ in powers of K we get 

2(3-2y) K 


He cx 


unp 


-c^ 

2 *^3,unp 


"0,unp 


-c^ 

2 *^2, unp 


)+o(a7£) . 


(97) 

y i-y \ 

If K is not too small, e.g., K 1/3 for ~ 1/10 and — A^ 3> —A T„ , as it is the case in 

the present fixed target experiments, we realize that cf u^p is in fact not numerically suppressed 
in comparison to cf ^,^p. So we would expect a contamination of the leading twist-two prediction 
by a A^/suppressed term, which, however, contains only twist-two CFFs. Even if this 
1/Q-contamination is small, Aq still fails to extracts solely the coefficient of the interference 
term and gets additive correction from cf ^^p, which stems from the gluon transversity (and higher, 
more or equal to four, twist effects). Since the gluon contribution is suppressed by a power of 
the strong coupling it may affect the twist-two coefficient ef unp in a modest way, though. 
Therefore, as a first order approximation the definition (|^) can be used for an order of magnitude 
estimate of the effects. 


^Note that experimental results from polarized deeply inelastic scattering give only a constraint for the for¬ 
ward limit of the antiquark-gluon-quark content of On the other hand, theoretical arguments that flavor 

small multi-parton contributions rely on model assumptions In our opinion, one cannot exclude that these 
correlations might be important. Thus, one should not rely on the WW-approximation and rather search for their 
signature in data. 
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The (definite charge) beam-spin asymmetry on an unpolarized target 


dlsL = 


^ , da^ — da^ /’^tt _ dfjl 


d(j) 


d(f) 


d(j)- 


dcf) 


,, da^ + da^ 
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(98) 


does not separate the interference term alone. So it does contain contributions from the squared 
DVCS amplitude: 

A ^ A2 _ 

(99) 
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Asl oc ± ^ /r,n(2^/^)s^,unp “ ^^^5° 
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with 


y„(2A7j) = -- 




sin(n0) 


(■ 22 ^ sin(n0) 
d(p- 




( 100 ) 


2/2 \Jo "^Pl(0)P2(0) ■ 

In the collinear limit, the asymmetry va ni shes. While for ?/ —> 0, it is determined by the twist-three 
coefficient The expansion with respect to K reads: 


2(3-22/) K 


^SLOCSi_„,,p 3^2-2/) 1-2/ 


(1 - 2/)(2 - 2/)-^^ „DVCS I \ /iQiN 

- ,ej . (101) 


The 1 / Q-suppressed effect stemming from the BH propagators induces a contamination by the 
second harmonic suppressed by K/{1 — y). Depending on the kinematics and the size of multi¬ 
particle contributions, this contamination together with may not allow clean access to 

the twist-two GPDs even from high-precision measurements of this asymmetry. Note that for 
general reasons, sf u^p is absent in the unpolarized interference (squared DVCS) term. 

The normalization of Asl is affected by 1/ Q-effects in the interference term, however, mainly due 
to the coefficient cjunp. 


5.2.1 Facilities with positively and negatively charged lepton beams 


As we have already mentioned, when lepton beams of both charges are available, this provides a 
clean separation of twist-two and -three GPDs. In these settings one can discuss charge-odd and 
-even parts of the cross section (^^, which extract the interference, and squared DVCS and BH 
amplitudes, respectively. The integrated charge-even part of the cross section does not contain any 
twist-three contributions — the interference term cancels there, while after azimuthal averaging 
the Co coefficient of the squared DVCS and all harmonics of squared BH amplitudes survive. So 
we will use it as a unique normalization of the asymmetries discussed below. Namely, independent 
from the target polarization we introduce 




^+^unp ^ ^-^unp 

d(j) 
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^BH^unp _j_ ^DVCS^unp 

dcj) 
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• Charge-odd part: 
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In this case, we end up with the interference term alone. However, because of the (/)-dependence 
of the BH propagators we have to include an additional weight factor and use the measure 

Pl(0)iP2(0)d0 


dw = 271- 


(103) 


SS’W)r,{4>V4‘'' 

for the azimuthal integration in order to compensate for the strong (/)-dependence of the product 
of lepton propagators. The measure dw has the properties 


r27r 

/ dw = 271, dw {(j)) = dw{—(j)) = dw{(j)27r). 

Jo 


Now we can exactly separate the Fourier coefficients in Eq. (pTl) . Note that dw has its minimum 
at 0 = TT, when the outgoing photon lies in the lepton scattering plane. In the case when its 
momentum is collinear to the lepton beam, this minimum approaches zero in the massless limit. 

To project out different harmonics, on can either (i) do the azimuthal averaging with appro¬ 
priate weights, namely, 

cos(n0)drc, sm{n(j))dw , (104) 

where n = 0,..., 3, or, (ii) use the fact that there is only a hnite number of terms in the Fourier 
sum of the cross section and integrate over different partitions of the azimuthal sphere. For n > 0 
this has an advantage of having a numerical enhancement by the factor of 4/7r with respect to 
the first method. Let us present the charge odd asymmetries (CoA), which distinguish the cosine, 
cos(n0), and sine, sin(n0), harmonics. 

The cos-harmonics, i.e., coefficients, are projected out by means of the integrals 
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(105) 

+ 2C!oAf(3), (106) 

(107) 


k=l 


/{2fc-3)7r/(2n) d(j) 

with n running over n = 2, 3 in the last equation. Here we adopt the analogous decomposition of 
asymmetries as for the cross section in Eq. (|3^ ) so that A = {unp,LP,TP}. The projection can 
be achieved by an appropriate flip of the target polarization vector. For a transversely polarized 
target, asymmetries are given in terms of two different combinations of CFFs. They are separable 
by the projection of the first odd and even harmonics in (p, while the average dip CoA"*"^ = 0 
vanishes. 

Next, the sin-harmonics, s^, can analogously be separated with the help of the formulae 


coAy, = A/; 
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i-k-K/n _ rl~(7^ 
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k=l 


' (fc—l)7r/n 


dcj) 


where n = 2,3. 


Charge-even part: 


Furthermore, we dehne the azimuthal asymmetries of the charge-even part. To do this in the 
cleanest way, we subtract at hrst the BH term. This might be possible in practice in the analysis 
of experimental data since the BH cross section is known exactly (up to electromagnetic radiative 
corrections), see section E3. with the nucleon form factors measured elsewhere. This gives us the 
squared DVCS amplitude 

2 ^ dvcs ^ = d~^a -|- d~a — 2d^^a . 


Then an appropriate azimuthal averaging, now with the conventional measure dcj), separates the 
cos-harmonics, via 
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with n = 1, 2, and sin-dependent coefficients by means of 
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with n = 1, 2. 

To conclude, an experimental facility having electron and positron beams is an ideal place to 
study GPDs and, thus, to extract the fundamental information about the spin structure of the 
nucleon. 


5.2.2 Facilities with single-charge lepton beam 

Now let us come to the situation when only one kind of the lepton beam is available. Here the 
study of single (lepton or nucleon) spin asymmetries allows to remove the background BH cross 
section. Note, however, that when both the beam and target are polarized, and one studies double¬ 
spin asymmetries, one gets the contamination from the BH harmonics too. In the single (lepton or 
hadron) spin experiments, still the twist-two coefficient sf is contaminated by power-suppressed 
effects, since both the interference and squared DVCS terms contribute. The best one can do in 
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these circumstances is to cancel completely the twist-three part of the interference term in the 
numerator. However, still one will have the presence of the power-suppressed DVCS cross section. 
For instance, for the single lepton-spin experiment one enables to probe plus (1 — y)A‘^/yQ‘^ 
corrections from |7 dvcsP and this can be done with the formula 


SSAi = 


, da^ — da^ r'^'^ , da^ — da^ 


dw 


d(j) 


dw- 


d(j) 


(113) 


iB-l) 


k-\-l 


rkir/6 _ ([(jl 


k=l 


' (fc— 1 ) 71/6 


dw- 


dcj) 


/■2’r da^ -|- da^ 

L 


Analogous extraction of the twist-two coefficient sf Lp is available for the nucleon-spin asymmetry 
(with the unpolarized lepton beam). The projection of the same components can be achieved by 
weighting the integral with sm{(j))dw. 

For smaller value of y the contamination from the squared DVCS term may be large. Let us 
demonstrate that a separation of the interference and squared DVCS terms can be achieved by a 
Fourier transform. The multiplication of the cross section with dw/dcf) induces new harmonics in 
the squared DVCS term. Projection of all three harmonics, i.e., measuring also 


2n 


SSA„ = ( i:(-l)‘+' r'" 

\k=l J(k-l)n/n dcp 

(where n = 2, 3) provides the desired Fourier coefficients: 
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OC SSAg, (115) 


with Ws^km = dw sin(fc0) sin(m0). Thus the extraction of twist-three harmonics is, in principle, 
possible, however, it requires high precision data. 

Such a modihed Fourier analysis might be employed to separate the coefficients of the inter¬ 
ference and squared DVCS term also in unpolarized and double spin-flip experiments. One has to 
be aware that the weighted cross section contains now four odd and hve even harmonics and that 
the BH cross section must be removed. 


6 Properties of Compton form factors 


To compare twist-two and -three contributions quantitatively, we use essentially the same ansatz 
for the GPDs as in . GPDs are constrained by sum rules and reduction formulae to the forward 
kinematics, which give us a guideline on how to model them. One may hope that these constraints, 
once implemented into a parametrization, provide a realistic order of magnitude estimate for cross 
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sections and asymmetries. However, we shonld emphasize that onr evaluations (as well as by 
other authors) are strongly affected by the model ambiguities involved. Experimental data will 
necessarily constrain GPDs via theoretical formulae. Therefore, it is important to understand the 
(i) theoretical uncertainties and (ii) influence of the model parameters on predictions. While the 
first issue is relatively straightforward to handle, provided one has enough knowledge on higher- 
order and -twist corrections, the second issue relies, more or less, on the experience that has to 
be collected by applying different GPD models. 

6.1 Ansatz for GPDs 

To give a quantitative estimate of the GFFs we have to model, or parametrize, GPDs. Let us 
discuss this issue for the twist-two sector, where we set rj = —Note that the time-reversal 
invariance and hermiticity imply that the twist-two GPDs are real valued functions symmetric in 

V- 

6.1.1 General remarks 

For the reader’s convenience let us remind the definition of the twist-two GPDs in terms of light-ray 
operators and their reduction to the conventional parton densities. We will use this information 
for the modeling of the former. The twist-two quark GPDs are given by 

{P 2 \' 4 ^{—Kn)'y ■ n'ip{Kn)\Pi) = J ■ n A"^) + e ■ nE{x,^, > (116) 

{P 2 \tl^{—Kn)'y-n'y 5 'il:{K,n)\Pi) = J ■ n H{x,^, A'^) + e ■ n E{x,^, A'^)'^ > (117) 

for the vector and axial channels, respectively. The Dirac bilinears are introduced in Eq. (^. We 
have contracted the free Lorentz index with the light-like vector in order to render the operators 
twist-two. For brevity we dropped the path ordered exponential between the quark fields, which 
makes the nonlocal operators gauge invariant, and did not indicate their scale dependence. This 
dependence arise from the renormalization procedure of the light-ray operators and is governed 
by a renormalization group equation Q. 

On the other hand unpolarized and polarized quark (/, A/) and antiquark (/, A/) densities 
are defined similarly by equations 

(p|' 0 (O )7 ■ ?7,'0 (k?7,)|p) = 2p ■ n J dx , (US) 

{p\2p{0)-f ■ nj 5 'ip{Kn)\p) =2s-n dx {Af{x)e-^^^^P-^^ + Af(x)P^^^P-^A , (119) 


29 


with the nucleon polarization vector = |f/(p) 7 ^ 75 f/(p). Here, the limiting procedure from the 
off-forward kinematics is dehned as A ^ 0 , so that Pi = P 2 = p, and one gets a restriction on the 
GPDs H and H: 


H{x, 0,0) = f{x)9{x) — f{—x)9{—x ), H{x, 0,0) = Af{x)9{x) + Af{—x)9{—x). (120) 

The hrst moments of the twist-two GPDs are equal to the corresponding parton form factors in 
the nucleon 

j'^dxH{x,^,A^) = Fi{A^), j'^dxE{x,^,A^) = F2{A^), ( 121 ) 

J^^dxH{x,^,A^) = Gi{A^), J^^dxE{x,^,A^) = Gp{A‘^), 

i.e., Dirac, Pauli, axial, and pseudoscalar form factors, respectively. 

For the gluonic GPDs we have the parametrization 

{P 2 \nA^{-Kn)G^ 0 {Kn)nfi\Pi) = ^ (P ■ n) £^dx 

X {h -nHc^x^i^A^) -1- e • nP g(x, A^)} , (122) 
{P2\naGapi{-i^n)iG^,0{Kn)np\Pi) = ^ (P ■ n) j 

X ^ - nHcix.i.A^) + e-nPG(a;,^, A^)} , (123) 

for the even- and odd-parity sectors, respectively. To make the contact with the forward parton 
densities, let us recall the definitions of the gluonic distributions, 

{p\naGa^,{^)G^^{^^n)n 0 \p) = (p ■ nf dx xg{x) {e"“''^+ -F , (124) 

{p\naGat,{O)iGf, 0 {Kn)n 0 \p) = {s ■ n) {p ■ n) f dx xAg{x) . (125) 

Therefore, we must have 

Hg{x, 0, 0) = xg{x)9{x) —xg{—x)9{—x) , Hg{x, 0, 0) = xAg{x)9{x)+xAg{—x)9{—x). (126) 


The first moment of these GPDs is related to gluonic form factors. The vector form factor is 
measurable in diffractive meson production and was theoretically estimated in Ref. . 


In different regions of the phase space GPDs share common properties with conventional 
inclusive parton densities for |a:| > and exclusive distribution amplitudes for |x| < Moreover, 
their experimental exploration will resolve the target from different perspectives and will enable 
to obtain its real-space image ^ . 
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To model the GPDs we will use the representation in terms of the so-called double distri¬ 
butions (DDs) m, In this representation one treats both s- and t-channel momentum flows 
independently with the corresponding fractions y and z, respectively. For the matrix element of 
the composite operator constructed out of scalar field operators, (/>, we would have 

{P 2 \ 4 >{-Kn)ct>{Kn)\Pi) = dy /{y, z, A^) 

7-1 7-l+|j;| 

= (127) 

The function introduced on the first line of Eq. (|127|) is the aforementioned double distribution 
f{y,z,A‘^), and g(x,^,A^) is the corresponding GPD. The expression of g(x,^,A^) in terms of 
f{y, z, A^) is given via 


q{x,^,A‘^)=[ dy ( ^ ^ dz5{y P^z-x)f{y,z,A^), 

7-1 7-l+ly| 


(128) 


and the inversion is known as the Radon transformation |^, |5^ 


However, in the case of spin-1/2 quarks there is an additional vector index associated with the 
Dirac matrix inserted in between the field operators. This leads to the appearance of extra Lorentz 
structures in the decomposition of the matrix element, which once discarded, lead to a violation 
of the polynomiality condition for the GPDs. In the historically hrst dehnition, the relevant 
structure was erroneously neglected and one used the same Eq. (p.28|) even for the nucleon matrix 


elements of QGD operators. This problem has been resolved by taking into account an extra term, 
the so-called D-term [Q. It is only concentrated in the exclusive region, i.e., \x/^\ < 1. Such 
contributions arise from isolated mesonic-like states^. However, there is yet another possibility to 


resolve the polynomiality problem |^. Unfortunately, in this case one has to deal with a more 
singular double distribution at j/ = 0, which requires, in general, a regularization that is invisible 
in the forward case. We will accept for the time being the former solution, which allows us to 
have a comparison with numerical results of Ref. []^ . 


6.1.2 Models 

As we have said above, we assume a two-component form for the GPDs. Namely, for the unpo¬ 
larized case one has for each Ttype quark 

H,{x, e, A^) = g,(x, e, A^) + ^0(1 - |x/e|)D(x/e, A^), (129) 

_ 

^Frorn the partonic interpretation of GPDs, it is clear that the exclusive region is always governed by mesonic- 
like states. In general, there is a cross talk between the exclusive and inclusive regions, e.g., this ensures the 
polynomiality condition for moments |Q. If this is not the case, we call them isolated mesonic-like states. 
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( 130 ) 


Ei(x,(,A^) = A^) - -rp«(l - |i/f|)C(i/f, A^) 

iV f 


with the D-term contributing equally to all Nf active quark species. Similar equations hold for the 
gluonic GPDs with i = G and a gluonic D-term, Nj^D{x/^) —> Dg{x/^). Note that the D-term 
drops out in the sum H + E and, thus, has to be an antisymmetric function of x/^. For three 
active quark flavors we set Nf = 3. In the following we take an oversimplihed factorized ansatz of 
the A^-dependence from the other two scaling variables x and ^ for all functions. Note that due 


to the antisymmetry of the D-term it does not enter into the sum rule (|121|) and, therefore, its 
A^-dependence is not constrained by it. As a model we take 


D(z,&^)=(i-N 


-3 




D{x/i ), 


(131) 


where the A^-dependence is characterized by the cutoff mass tud, considered as a free parameter. 

Obviously, g* is a sum of valence- and sea-quark contributions. According to this we make a 
decomposition in these components and extract the momentum transfer dependence into valence 
and sea form factors, namely. 


A^) 


F i.val 
1 


(A2)j™>(i,J) + Ff‘(A2)C“(i.?). 


(132) 


for i = u,d, s, and = 0 for the s-quarks. 

At this point let us say that the (a:, ^)-dependence of the function r(x, A^) is not constrained 
at all. However, since its hrst moment is given by the Pauli form factor, we extract it from 
r(x, A^) and naively set the remainder of r equal to g(x, ^), i.e., 

r,(x, e, A^) = gr'(x, o + o • ( 133 ) 


The valence u- and d-quark form factors can be extracted from the proton and neutron form 
factors via the formulae 


r\ T7Tli,val 
^-^1,2 


(A") = + F”2(A") , 



(A") = FP(A") + 2F”2(A") . 


(134) 


The latter are known fairly well from experimental measurements and can be parametrized by 
dipole formulae in the small-A^ region 
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\ K, 


^ ll - N 


-2 
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G"(A2)=0. (135) 




Note that we have set the neutron electric form factor equal to zero since at small momentum 
transfer it vanishes as a hrst power of A^ with —1/6 times the neutron charge radius as a 
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coefficient, which is quite small, ~ —0.113 fm^ . Here the proton and neutron magnetic moments 
are K,p = 1.793 and = —1.913, respectively. We have used Sachs electric and magnetic form 
factors related to Dirac and Pauli ones by 


. (136) 

G‘„(A") = F;(A") + F‘(A=), (137) 


respectively. They are characterized by the cutoff mass my = 0.84 GeV, see, e.g., [^. The 


sea-quark form factors are guided by the counting rules and read in terms of Sachs form factors 

A2 


G sea 
E 
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1 + ^Si 


^sea 


(A^) = 1 


mf 


(138) 


with yet another mass cutoff msea- The slope of is given by Bgea = 5F®®®'/5A^|a2=o = 
3/mgga ~ Ksea/4M^. The two free parameters Bgea and Ksea will be specified below. Note that Ksea 
enters in the sum rule that gives the fraction of the orbital angular momentum carried by quarks 
0. With our ad hoc assumption r* = g* this sum rule reads 

•I’-iim E \l'+ (139) 

i=u,d,s ^ 

= — |(1 -|- Kp -|- Knl‘F) -|- (1 -|- Kp -|- 2/D„) + (1 + Ksea) P®®®'! , 


where the momentum fractions P* carried by quarks can be deduced from deeply inelastic data 
alone. 

The D-term does not affect the forward limit of GPDs. Therefore, we use Eq. (|128|) and model 
the quark fi{y, z) according to Ref. as a product of a prohle function vr with the conventional 
type-i quark fi{y) and antiquark fi densities. Namely, for valence and sea contributions we have 


fr\y^ = fr\y)^iy)^i\y\^ F Kai ), (i40) 

z) = {fiiyWiy) - fi{-yW{-y)}7r{\y\,z-, b^ea ), (i4i) 


with the usual definition for the valence density /™*(2/) = fi{y) — fi{y)- Here the profile is assumed 
to be universal for all valence- and sea-quark species and reads 


'wiy.z^h) 


^ ('’ + !) |(i-»)^-^t 

yFr(6 + i) (i-!/)^^ 


(142) 


This ansatz ensures that the reduction formula (|120|) holds true. The parameter b encodes the 
skewedness effect, i.e., a larger value of b suppresses the .^-dependence, and g*(x, reduces 

to the parton density f^{x, Q^) in the limit b —>■ oo. This limit will be called the forward parton 
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distribution (FPD) model |Q. The idea driving the form of the ansatz (|142|) is that the profile 
function 7r®(|/ = 0 ,z) mimics a mesonic-like two-parton state characterized by the longitudinal 
momentum fraction In the case of valence quarks one may choose the asymptotic distribution 
amplitude, which implies b^ai = 1. In the case of sea quarks we consider bsea as a free parameter, 
which will be discussed below. It will be adjusted by the comparison of our predictions to the HI 
small-XB data. 

We similarly assume a factorized ansatz for the gluonic DD 


where 


fo(y,z,A^) = Fc(A^)fa(y,^), 


fa{y, 2 ) = y{3{y)<)(y) - g(-y)<)(-y)}Tr(\ylz;bc) 


(143) 


(144) 


The quark D-term was given in []^, making use of the results of the chiral quark soliton 
model 1^ at the very low normalization point /i^ ~ 0.36 GeV^, in the form of an expansion in 
terms of Gegenbauer polynomials, — The expansion 

coefficients have been fixed under the assumption that the polynomiality can only be restored 


with the help of the D-term and fit to the chiral soliton model evaluation. Gompared to Ref. 
we can adjust its fall-off with in Eq. ([I31|) by the cutoff mass tud. Note that the evolution 
upwards to our input scale Qq = A GeV^ from 0.36 GeV^ reduces the absolute value of the 
Gegenbauer polynomial coefficients by about 25% or even more, provided that the gluonic D-term 
is neglected. This reduction is larger (smaller) for positive (negative) expansion coefficients of 
the gluonic D-term. However, to study the significance of such a term for DVGS observables, we 
ignore any evolution effects. 

The twist-two GFF 7Y(^, A^) at LO in the QGD coupling constant reads in terms of the model 
(|128|) with the D-term included, i.e., Eq. (|129| ), 
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2 Fr(A") {M:y)0(y) - M-y)0(-y)) ir(|9|,;; M 


where Nf = 3. Obviously, the replacements D{x/^) —> —D{x/^) and Fi F 2 lead to the 
definition of S. 

In the following we take the MRS A’ [Q parametrizations at Qq = A GeV^. This input scale 
ensures that evolution effects are small and are, therefore, not considered. We will neglect SU(2) 
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breaking in the sea and its charm component, i.e., u = d = s/2. We use these forward parton 
densities in our LO analysis of DVCS and we are aware that this may result in a deviation of 
our predictions from the forward structure functions, measured experimentally in deeply inelastic 
scattering. 


Xb 


F2ixB, ~ 2xbFi{xb, = XB,r] = 0, = 0, Q' 


TT 


(146) 


up to 20%. Here we employed the Callan-Gross relation, valid in LO of perturbation theory. 
Recall that for a general two-photon process in the light-cone dominated region with a virtual 
outgoing 7 -quantum, the OFF is given by A^, Q^). We will demonstrate below that the 

error involved in this approximation is small in comparison to the uncertainties induced by other 
sources. 


Note that the fraction of the orbital angular momentum, given by the sum rule (1^) with the 
assumed MRS parametrization, results into 


^ 0.3 + O.OOKsea , for = 4 GeV^ 


(147) 


Estimates based on lattice calculations and QCD sum rules provide J'^ = 0.3 ± 0.14 and 
= 0.15 ± 0.13 at /i^ ~ 1 GeV^, respectively. So we conclude that for consistency of these 
predictions, Kgea may vary in a wide range —3 < Kgea < 2. 

By the same token as above, the polarized quark GPD H is modeled according to 

H,{x, e, A2) = Gr\A^) Aqr\x, 0 + Gr(A2) Aqrix, 0 , (148) 

for i = u,d, s, and Ag™^(x, 0 = 0 1°^ 1^® s-quarks. Note that there is no D-term in this function 


17|. The polarized valence- and sea-quark contributions are 


Afr\v,z) = Afr\v)0{vMv\,Zib „{), 

^/r'iv.z) = {Aft{y)e{y) + Aft{-y)0(-y)'^ n{\y\, z;b^,) , 


(149) 

(150) 


with the usual definition for the polarized valence-quark density AfA^{y) = Afi{y) — Afi{y). 
Obviously, A/J^* = 0. In the forward limit this formula leads to Eq. (|120|) . For the forward 
densities we use the GS A parametrization with an SU(3) symmetric sea, again at the input 
scale = 4 GeV^. 

The axial form factors are 

-2 , „ , _3 


Gr'(A2)= 

\ '^Ay 


GfOAO = (1 - —A 


(151) 
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where i runs over u- and d-quark species and the scale is rriA = 0.9 GeV H- The slope Ba is 
considered as a free parameter, which should be hxed from experimental data. 

The polarized CFF A^) in terms of the model ( 128 ) has the form 
~ - 1 B dx rl rl-l 


= 7 


/ dy dz6{y + ^z — x) 

J-l 


(152) 


{ i-i 1 - x/( - iO J-i i-i+lsl 

X E 0dGr'(A=)(A/--(!,)9(!,) + Ar>(-!,)9(-;,))^(|!,7;6,.,) 


i=u,d,s 


+2Gr(A=) {Mi(y)e{y) + A/,(-9)9(-j,)) 2 ; 6. 


For E we accept the pion pole-dominated ansatz of Ref. 

A^) = -EHxAA'‘) = 

with F^(A^) taken in our estimates in the form given in Eq. (39) of Ref. [p^, 


(153) 


F.(A2 


/ 1.7(A^-m^)/GeV^ ^ 

m2 - A2 (1-A2/2GeV2)2 


(154) 


valid for |A2| M'^. The pion mass is m^r ~ 0.14 GeV. For the pion distribution amplitude 

we take for simplicity its asymptotic form 4’t,{u) = =6u(l — u), 0 <u < 1. 

Finally, the parity-odd gluonic GPD, Hq-, is modeled via the reduction to DD, and the latter 
taken in the form 


A/g(|/, z) = y {^g{y)e{y) + ^g{-y)e{-y)} n{\y\,z- he) ■ (155) 


6.2 Features of Compton form factors 

Now, before we come in section to numerical estimates for the cross sections and asymmetries, 
it is worthwhile to study the properties of GFFs in the twist-two and -three approximations: the 
magnitude and shape of GFFs, the ratio of their real to imaginary part, the relation between 
twist-two and -three GFFs in the WW-approximation, and the effect of NLO corrections. Within 
our approximation it is consistent to take ^ = x-q/{2 — xb)- Since we use a factorized form of 
GPDs with A2-dependence separated from the scaling variables, we can safely set A2 = 0 in the 
most of estimates made in this section. These calculations are done with = Qq = 4 GeV^. 


6.2.1 Twist-three versus twist-two effects 


As observed above, the real and imaginary parts of the partonic GFFs are convex or concave 
functions of ^ (see the archive version of |P3). In Fig. □ we show the GFFs for the valence 
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up-quarks (a,b) and up-antiquarks (c,d) based on the MRS A’ parametrization without the D- 
term and with 6vai = Ken = 1- Here we rescaled the CFFs with and used = 0.441 and 
Q,sea _ which determine the small momentum fraction behavior of the corresponding forward 
parton densities [Q. The D-term affects only the real part of the CFFs by a shift, e.g., and 
are shifted by a constant value of ~ —1.66 and ~ —0.22, respectively. Then, the 
entire CFF = J2i=u,d,s'^ii.O have a signihcant change: it starts now with a negative 

value at large ^ and turns over to a positive one at xb ~ 0.2. Obviously, the sign of 3?e7-f in 
the valence-quark region]^ is determined by the cancellation between the positive sea-quark and 
negative valence-quark effects, both for double distribution ansatz alone and with the addition of 
the D-term. We will demonstrate below that a clear signature for the D-term can be washed away 
at 7 ^ 0 since the A^-dependence of the sea quarks and the D-term is not constrained by the 
sum rules (|121|) as we have discussed in section |6.1.2| after Eq. ( p.30|) . 

The shapes of the CFFs can be understood in the following way. The phases = Arg 
and the absolute values \T-Ci{^) \ are monotonically increasing functions of the scaling variable K H 
is interesting to note that in the WW-approximation the phase difference A0j = 0i(O 
is not very large. Neglecting the D-term and 7^^, we hnd that |A0™^| < 0.4 and < 0.16 for 

all values of K Moreover, these differences will vanish in the limit ^ 0. This property is obviously 

signihcant for observables that arise from the interference of twist-two and twist-three CFFs. Note 
also that the absolute value of the effective twist-three CFFs in the WW-approximation is of the 
same order as the twist-two ones. Thus, there is no numerical enhancement of the twist-three 
effects in this approximation. The relative size of twist-three contributions as compared to the 
twist-two harmonics is governed solely by kinematical factors appearing in the Fourier coefficients 
computed in section 0. 


6.2.2 Small-XB behavior 

Let us now study the properties of the CFFs for small xb in LO approximation. At intermediate 
momentum transfer this kinematics was discussed within the BFKL approach in Ref. II64II . In 
Fig. I we observe a scaling power-law behavior of the CFFs for small which has been already 
established in Refs. |65| , |66| and the archive version of for the FPD model. The small-.^ 
behavior of the CFFs is governed by the small momentum fraction asymptotics of the parton 
densities. If the latter behave like with a > 0, then the real and imaginary part of goes 
like and the ratio of the real to imaginary part is given by tan {{a — l)7r/2) = — cot {an/2). 
The small x-behavior of polarized parton densities, i.e., y~°‘, induces an analogous growth of H, 


®It implies that xb ~ 1/3, which we understand in the following as the region 0.1 < xb < 0.6 . 
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Figure 2: The real (a,c) and imaginary (b,d) part of the u-quark CFFs in LO (solid) and NLO 
(dash-dotted). It displays also the twist-three functions in the WW approximation ^“ 7 -^efr-ww 
with = 0 (dashed), and the additive contribution (g) H (dotted), which enter the CFF 

multiplied by a kinematical factor singular at The valence and sea u-quark 

contributions are shown in (a,b) and (c,d), respectively, as a function of xb with xb > 0.05. The 
GPDs are taken at A^ = 0 and = 4 GeV^ with 5vai = &sea = 1- 

however, the ratio of the real to imaginary part is now tan (57r/2). 

Such a behavior can be understood in a model-independent way. The essential assumption is 
that the double distributions, which belong to the class of mathematical distributions, develop an 
expansion in the vicinity of the point 1 / = 0 of the form 


fi{y,z) = y "‘{di( 0 ,z)+ ...}. 


(156) 


Note that a (5-hke singularity a.t y = 0 gives a contribution that is entirely concentrated in the 
exclusive region. We refer to it as an isolated mesonic-like state. In the forward limit of GPDs, 
such a contribution appears as a 5-hke singularity and, thus, it will affect sum rules for parton 
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densities. So we may conclude that such singularities are absent in the spin non-flip GPDsQ H and 
H, however, they may be present in the spin-flip ones. Furthermore, we assume that f{y,z) and, 
consequently, also d(0, z) vanish fast enough when they approach the boundary of their support. 
By ‘fast enough’ we mean that |(i(0,2;)| < |1 =F with £ > 0 for ^ ±1. We only 

discuss the region 0 < ?/ in the following. The results of antiquark contributions for negative y 
can be easily deduced from the sea quarks at ?/ > 0 employing the (anti)symmetry property of 
corresponding GPDs. We also assume that Oj < 2 in the vector and dj < 1 in the axial-vector 
sector to ensure the existence of certain integrals. 

Since GPDs are reduced to the parton densities fi{y) according to Eqs. ( |120|) , we see that at 
small momentum fraction y, the parton densities are given by j\ dz di{0, z) with the expo¬ 
nential ttj, which determines their small-y behavior. Furthermore, a straightforward calculation 
shows that the imaginary and real parts of the twist-two GFFs at LO behave like: 

Q^) = Q^r”‘ £ dz . (157) 

SeWi(5. A^ Q^) = tan ((a, - 1)0 SJmWi(5. Q"), 

where the parameter a* may depend on and Q^. The first equation agrees with the earlier 
considerations in Ref. and the second one was discussed in Ref. |^. An analogous formula 
holds true for the spin-flip contributions Note that the D-term is not important in this limit, 
however, the small-^ behavior can, in principle, be altered by other terms concentrated at y = 0, 
i.e., -^5{y). So far such a behavior is not excluded by sum rules. 

Let us also give the ratio of the imaginary part of the GFFs to parton densities, 

where = Q:i(A^ = 0). For the GPD ansatz, we are using with 6gea = 1, the leading be¬ 
havior comes from the sea quarks with Ogea = 1-17, which provides their ratio ( |158[ ) of order 
~ 1.757rQ^F®®'^(A^) for each sea quark species. 

For the axial-vector channel, analogous relations can be derived, however, with a difference 
that the ratio of the real to imaginary part is now tangent rather than minus cotangent, 

A^ S=) = dz , (159) 

3fJe7Yi(^, A£ = tan Am7Yj(^, A£ . 

^Unfortunately, this is not necessarily true since they can vanish in the forward limit with A^, too. 
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The ratio of the imaginary part of Tii to the polarized parton density Ag* is determined by 


Ag.(e, Q2) 



!\dz (ii(0,z,0, Q2) 


(160) 


with 5° = 5j(A^ = 0). However, in contrast to the vector channel, the small-xe behavior of 
the polarized parton distributions has not been explored in experiments. Usually, one assumes a 
similar behavior of the polarized valence- and sea-quark distributions with dj 0.5. Thus, the 
leading behavior of is determined by all parton species and in the forward limit it is suppressed 
with respect to by an extra or similar. 

In the case of twist-three contributions, the power behavior remains the same, however, the 
normalization of the real and imaginary part will be changed. It is sufficient to discuss the 
convolution with the -kernels. The behavior of the effective twist-three contribution in the 
WW-approximation of Eq. (|8^), can then be simply predicted. For the imaginary part we hnd a 
modihcation of the integrand by the hypergeometric function: 


3m 


(do? ® •f’-}= 

iXi 


-OLi 


/-I (1 — z)°‘i 


2-^1 


1, CXi 

1 -|- Q-i 


1 + ^ 
l-Z 


(161) 


The ratio of the real to imaginary part remains the same as in the twist-two case, i.e., — cot(7rQ;/2) 
and tan(7r5/2) for the convolution with and kernels [see Eqs. ( 157 ) and (|159|) ], respec¬ 
tively. 

From Fig. ^ we observe that the relative deviation of the dominant partonic CFFs and 
from their asymptotic values given in Eqs. ( p.57|) and ( p.59| ), respectively, is of order of 20% for 
xb = 10“^ and decreases to 10% level, and even further at xb < 10“^. Note, however, that a faster 
convergence of the effective twist-three CFFs to its asymptotic form (|161D in comparison to the 
twist-two sector is caused by an essential cancellation occurred in the WW-approximation. For 
instance, the convolution ® alone has a large deviation of order of 40% for xb = 10“^. 

Finally, let us add a word on the isolated mesonic-like states, which we mentioned at the 
beginning of this section. It is easy to show that these 5-like singularities, i.e., ^6{y), convoluted 
with the hard scattering part yield contributions proportional to Since the hard scattering 

part possesses dehnite symmetry properties, the exponent n is odd (even) for S (T). 


6.2.3 Radiative corrections 


Now we address the theoretical uncertainties due to radiative corrections, which have been ana¬ 
lytically worked out in a number of papers for the twist-two sector at NLO [0 0 0 0, 

and numerical studies have been performed in Refs. |^, ^ . For the quark sector alone the 
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Figure 3: The relative deviation of the sea u-quark CFF from its asymptotic value is displayed 
for ‘^mH and 3fJe7-f/9m7-f as the solid and dashed line, respectively, in (a). The same quantities 
are plotted for ^ g as dash-dotted and dotted lines, respectively. The right 

panel (b) shows the same as in (a) but for the valence u-quark CFF 7i. The GPD-parameters are 
the same as in Fig. 


NLO amplitudes are shown as dash-dotted curves in Fig. The one-loop corrections turn out to 
be moderate for each quark specie. However, model-dependent considerations that involve glu- 
onic contributions showed that NLO effects can be quite sizable both in the valence ^ and 
small-XB 1^ region. In the following we have a closer look to this issue and propose a solution to 
the problem of large NLO effects. 

The NLO corrections are evaluated with the same set of GPDs as earlier. We disregard the 
evolution of the coupling and set it equal to a^/vr = 0.1. The D-term and pion-pole contributions 
will separately be considered below. If the real part of a GFF is close to zero at some e.g., 
due to cancellations between different quark species, the relative radiative corrections are blowing 
up. Thus, it is more convenient to discuss their magnitude in terms of the absolute values |JF| = 
y(3fJeJF)2 -I- (3mjF)2 phases Arg JF. 

Let us study at first the naive scale setting condition fip = Q- In this case the radiative 
corrections to the absolute value of the GFFs, i.e.. 


^abs _ |7^NL0|/|^LU| _ 


LOi 


(162) 


are rather moderate provided gluon contributions are totally neglected. Then for the choice of 
the DD-model 6-parameter 6sea = 1, one observes that is of order —10% in both valence and 
small-a:B region, as shown in Fig. ^ (a) by the dashed line. For the FPD case when 6sea —^ oo, 
the is a factor of 2 larger as compared to i.e., ~ —20%. The radiative corrections 
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Figure 4: The relative NLO corrections IZ (see Eq. |162|) to the absolute value (a,c) and the phase 
(b,d) of the CFFs Ti (a,b) and Ti (c,d). The magnitude of radiative corrections with neglected 
gluon GPD contributions are displayed as dashed and solid curves for the factorization scale setting 
lip = Q and /ii? ~ respectively. The complete (quark + gluon) NLO result is given for 

these both scale settings by dotted and dash-dotted lines, respectively. The coupling is set to 
as/p = 0.1, the quark GPD-parameters are the same as in Fig. and be = 2. 


stemming from the quark sector alone to the phase of GFFs are defined as in Eq. (|162| ) but with 
the obvious replacements of by and Ti by AigTi. They are smaller than 20% in the 
valence quark region, and are on a percent level at small x-q, see the dashed line in Fig. ^ (b). 
Qualitatively, the same features hold for Ti. However, quantitatively we hnd here a slightly larger 
effect in the valence-quark region where relative corrections are of order of —30% and ±20% for 
the absolute value and phase, respectively, see the dashed curves in Fig. | (c) and (d). 

To have a feeling on the magnitude of GFFs we give a few numbers. For we have 

= 5 . lo-Q = 0, = 4 GeV^) ^ (0.39 ± 1.37Q ■ 10^, 
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which is considerably larger than its polarized connterpart, 


= 5 ■ 10-^ = 0, = 4 GeV^) ^ 200 + 163i, 


mainly reflecting the fact that 5* < Osea- 

So far, at the level of amplitndes the NLO radiative corrections in the qnark sector were 
moderate within the model with bsea. = 1- This gives us some confidence in the applicability of the 
perturbative approach to the DVCS process. Now we come to the gluonic contributions, which 
enter for the first time at order 0{as). For fj-p = Q they are displayed as dotted curves in Fig. |^. As 
we see, with our model for gluonic GPDs their effects in the amplitudes are large. However, the size 
of the product ag/^nx gluonic GFF, compared to the LO result, is not a measure for the quality 
of perturbation theory, it rather gives the ratio of ‘leading’ gluonic to quark GFFs. Of course, 
one cannot exclude that these large Os-contributions are merely generated by an unrealistically 
oversized ansatz for the gluonic GPDs. 

To understand the physical meaning of these large contributions, it is instructive to study the 
evolution of the gluonic GPD. We use throughout 6^ = 2 in the profile function (|142D. Taking 


the GRV parametrization at its intrinsic low scale Qq = 0.49 GeV^, it has been observed that the 
gluon GPD changes under evolution much stronger than their non-singlet quark counterparts 


This is of course also true in the forward case. However, if we take the MRS A’ parametrization at 
the input scale Ql = A GeV^, we observe a strong scale dependence of the NLO result, which is in 


line of observations made in Refs. |^, |^. Since the NLO coefficient functions can be expressed 
in terms of evolution kernels |^, thus, a strong change under evolution implies a large gluonic 
contribution in the one-loop approximation. Analogous to deeply inelastic scattering one can 
study the scale dependence by taking the logarithmic derivative of the GFFs, i.e.. 




:{Q^)Ei=uAsCl 


(-) 


1/59 (g) iP. + 


H, 


G 


K) 


WK) 


(163) 


where 1/99 and 1 / 9 G ^re the LO evolution kernels. This quantity is demonstrated in Fig. for two 
quite different factorization scale settings at xb < 0.1, where the gluonic term produces, indeed, 
a dominant contribution. 

Of course, it would be delightful if gluon contributions would be also moderate. Then DVGS 
observables can be estimated relying solely on LO approximation, and the phenomenology of 
DVGS can in the first stage be restricted to the quark sector alone. However, even in the opposite 
situation we can minimize the size of large radiative corrections, originated by strong evolution 
effects. The restriction to the quark sector can be achieved by an appropriate factorization scale 
setting, namely, in such a way that the gluonic contribution does not enter the amplitude, e.g., at 
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Figure 5: The real (a) and imaginary (b) part of for MRS A’ parametrization with 

&vai = bsea. = 1 at = 0. The dash-dotted (dotted) line shows the scale violation due to the 
evolution of quarks, while the solid (dashed) one includes also the evolution of gluons for the scale 
setting /ii? = 3.8Q (/i^ = Q)- The parameters are set as in Fig. 

NLO 

(y _ /*1 

^g(^, Q^) = ^ I] Ql J_ dxCG{x/^,Qyfil)HG{x,^,fil)i,^^i,uvcs = 0. (164) 

i=u,d,s ^ 

This scale setting procedure will go under the name of the DVCS scheme. It requires two {x, ^)- 
dependent factorization scales inside the convolution, i.e., one for parity even and odd CFFs, 
respectively. Consequently, this procedure provides two (one for real and one for imaginary part) 
real-valued effective scales that depend on for each CFF. 

An inspection of the imaginary part of the NLO coefficient in the vector case tells us that 
the desired factorization scale is almost ^-independent. For our model it is given by ~ 

\/2eQ ~ 3.8Q. Since the real part is expressed in term of the imaginary part by means of a 
dispersion relation, we expect that also the absolute value of Tic is minimized in this way. This 
is demonstrated in Fig. where we display the relative NLO corrections, see Eq. (|162|) for the 
dehnition, coming from quarks (dashed, solid) and from all partons, including gluons, (dotted, 
dash-dotted). We also employ the scale found for H in the evaluation of H in the valence quark 
region. Here we should mention that with our models, the CFF Ti can be safely neglected in 
and for an unpolarized target. In the vector (a,b) and axial-vector (c,d) cases the 
gluonic contributions drastically reduce the LO quark contributions for the naive /i = Q scale 
setting (dotted line). For = 3.8Q, the radiative gluonic NLO corrections are small (cf. dash- 
dotted and solid curves). The displayed change of the quark radiative corrections (for the absolute 
value from —10% to —20%) should be considered as a bound, since we did not include evolution 
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effects. Doing so, they will compensate the change we have displayed. Let us emphasize again, 
that our quantitative considerations are model-dependent. Taking, for instance, the FPD model, 
at = Q we observe the value of is twice larger in both sectors. For the choice /i^ = 3.8Q, 
the gluonic NLO contribution is now about 10%, while the quark and the whole NLO corrections 
are even smaller as demonstrated in Fig. ^ (a). 

Let us now come to contributions that show up only in the exclusive region. For the D-term 
we hnd to NLO accuracy 

^D-term = -4.98 |l + ^ 1^1.83 - 2.34 In + ■ ■-j . (165) 

If we set the factorization scale = Q, the relative correction is positive and of order of 10%. If 
the gluon model is realistic, we should better employ the DVCS scheme, in which the D-term is 
reduced by about 20%. 

For the pion-pole contribution the scale dependence drops out in the hard scattering part since 
we have chosen the asymptotic form of the pion distribution amplitude. We get at NLO 

= ( 166 ) 

which acquires the correction of order 15%. 

6.3 Summary 

Let us summarize the results of this section. We introduced a simple model with several free 
parameters based on the factorization of and (x, dependence. The latter is modeled 

according to the DD ansatz with MRS A’ and GS A parametrization for the parton densities 
at the input scale Qq — 4 GeV^. This oversimplification is only justified by the fact that there 
is no model-independent knowledge available about the functional dependence of GPDs on their 
arguments. 

In this model the NLO corrections to the partonic quark GFFs are moderate. The theoretical 
uncertainties entering at NLO mainly arise from the lack of information on the form of gluonic 
GPDs. Also in the case when a realistic gluon GPD provides NLO corrections of the same order 
as the quark contributions at LO (in we observed a different overall sign with respect to 
the unpolarized sea-quark contribution within our model), we can exclude them from the NLO 
analysis by choosing the DVGS scheme. Although this favors a rather large factorization scale, 
the NLO corrections in the quark sector remain moderate. However, note that the results for the 
FPD model of is a factor of 2 larger than the corresponding estimate with 6gea = 1- Therefore, 
NLO effects can be quite sizable, especially, in the squared DVGS term. In this product of DVGS 
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amplitudes they are twice larger than in the interference term, which is linear in the them. For 
a rough estimate of observables, which are affected in an uncontrollable manner by the models 
accepted for quantitative estimates, it is enough to work in the quark sector at LO. However, for 
extraction of model parameters from experimental data with high accuracy, a NLO analysis will 
be inevitable. 

Moreover, we found from general assumptions, that the parametrization of the CFFs T = 
J2i=u,d,s ill tiiG small-a:B region is rather unique, and is in fact governed by the small ^/-behavior 
of the DDs. The skewedness effect is included in the normalization of the CFFs. Our analysis 
suggests that we can take the following general parametrization, derived in the DVCS-scheme at 
a given input scale Qq > IGeV^, at small x-q\ 


3?e 
55m 
■ 3fJe 
55m 
JJe ' 
Q'm 


3fJe 

55m 


VK.A") 

WK.A^) 

VK.A^) 


I - cot (a ir/2) | ^ 

I - cot (/3 tr/2) I ^ 1 1 1 _ 

|tan(aV2)|^^,^2)j_S,^.,_ 

I tan (pr/ 2 ) j ^ | + | J} Mj(A^)r^'-‘. 


where we have a(A^ = 0) ~ 1.2 and a(A^ = 0) ~ 0.5 and p, p > 1 are positive integers that reflect 
the appearance of isolated mesonic-like states. So they provide the dominant contributions in the 
spin-flip CFFs and could overwhelm the small-XB behavior of the spin non-flip CFFs. Note also 
that because of the small- and large-^ behavior, and the simple shape of the (partonic) CFFs, it 
is obvious that one can extend the parametrization in the whole kinematical region by allowing 
a weak ^-dependence in the A^-factors and additionally in the phase, while the M-factors remain 
^-independent. 


7 Quantitative estimates for DVCS observables 

In this section we give numerical estimates for DVCS observables in the twist-two and -three 
approximation at LO of perturbation theory for several choices of kinematical invariants and 
compare them, whenever possible, with available experimental data. To this end, we use exact 
expressions for the BH propagators and iC-factor, when the latter appears as a prefactor in the 
Fourier coefficients. In section o we use the DVCS data by HI collaboration to adjust free 
parameters in the sea-quark sector of our model. Then we consider asymmetries for unpolarized 
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fixed target experiments in section |7.2| . We compare our model expectations with the beam- 
spin asymmetry measurements of HERMES and CLAS collaborations and give then predictions 
for the E-91-023 experiment at Jefferson Lab IM. We deliver a general discussion about the 


magnitude of twist-three versus twist-two effects in these kinematical situations. Furthermore, we 
present a detailed numerical analysis of the size of twist-three effects for an unpolarized target 
and derive constraints for certain charge azimuthal asymmetries. Finally, we give in section |7^ a 
few estimates for longitudinally polarized fixed target experiments^. 

The reference frame used in the analyses by these experimental collaborations differs from the 
one we employ in our considerations, see Fig. by the direction of the z-axis. Both, HERMES 
and CLAS use the third direction pointing along the virtual photon’s three-momentum. Thus, 
their frame (distinguished below by the prime on symbols) is deduced from ours by the rotation 
around the x-axis on the angle vr. Therefore, the proton azimuthal angles in our frame cj) and 
their frame are related by 0 -f = 27r. The same relation holds for the real photon angles 
+ = 271. The proton and photon angles are related to each other by the equations = (j) + n 

and -|- tt in our and their frames, respectively. Therefore, the proton and real photon 

azimuth is given in terms of the angle 0 by 


07V — 271 — (f> , 


= 7T 


7.1 Small -xb estimates 


In this paragraph we address two questions: (i) What can we learn from DVCS in the small- 
xb region? (ii) Is there a kinematical window left to measure twist-three effects? Thanks to 
the small-XB parametrization devised in Eq. (|167|) and the analytical expressions for the Fourier 
coefficients, computed in section both of these questions can be answered in a straightforward 
manner. Although we only consider an unpolarized proton, the analysis is rather complex and we 
cannot cover all possible scenarios, which are allowed by the residual degrees of freedom left in 


the parametrizations (|167|) . 

Let us address the hrst question. We will now see how one can test the structure of CFFs 
and measure its small-xe behavior in a model-independent manner. This includes a test of the 


absence of isolated mesonic-like states, the measurement of Ti, and, perhaps, access to the spin- 

®We refer to the longitudinal polarization in our frame. In the laboratory frame with the z-direction defined 
with respect to the lepton beam it implies that the target has a transverse polarization either. The reverse is also 
true, if in the lab-frame the target is longitudinally polarized, in the frame Fig. we will have both longitudinal 
as well as transverse components of the nucleon polarization vector. The projection on the first one is easily done 
by integrating over the whole range of the angle y. 
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flip CFF £. The strategy is to rely first on the hypothesis of the sea-quark dominance in Ti and 
then to derive relations between different observables. A violation observed in such constraints 
will then be regarded as a signature of other contributions involved and not accounted for in the 
approximation adopted here. 

Employing Eqs. ( |167D , we can now drastically simplify the twist-two Fourier coefficients. For 
an unpolarized target the squared DVCS term [see Eqs. (^31) and (|69|) ] reads 

( 1 Nl / 1 


c, 


pvcs 

0,unp 


= 2{2-2y + y^ 


sm‘ 


{a 71/2) ^2“ 4M2 sin^ (/? 7r/2) 



(168) 


|^=xb/2 


where the term F^/A arises from the pion-pole contribution^ in S. The latter is absent in the 
interference term, which is given by [see Eqs. (^) and 




unp 
^l,unp , 


= 8K 


X 


Fi 


-i2-2y + y^)) 

M‘2-y) I 

— cot {a 7r/2) 1 

1 j ^ 


(169) 


A2 

4M2' 


— cot (/? 7r/2) 
1 



- IC=*b/2 

Note that the D-term can be safely neglected as emphasized above and Ti does not contribute 
if we assume that 5 -C a. The last assumption is expected to hold for parton densities in view 
of current experimental data. The normalization factors can be expressed in terms of model- 
dependent parameters. In our case they are determined by the partonic form factors, skewedness 
effects, and normalization of parton densities. Relying on the sea-quark dominance we hnd 


1V«(A= 


i=u,d,s 


2^sea(l “1“ T^sea) Tpsea./' \2\ -|- &sea *^sea)r(2 -|- 26gea) 




Qlo 


(170) 


:seaF(l -f- 6gga)F(2 -f- 26sj 

and analogous equation for with being replaced by Here 7sea contains the SU(3) 

flavor-breaking effect and for the MRS parametrization it is 7sea = 1/10. Note that the SU(2) 
flavor breaking, i.e., the non-zero difference d — u behaves as a valence-quark density and, thus, 
can be neglected in the small-XB region. In the MRS A’ parametrization the normalization factor 
for the sea-quark contribution is Ugea = 0.956 and a = ages. = 1.17. We should emphasize that at 
LO of perturbation theory, 

^ hm ^Nn{A^ = 0)= ^ g 2 ^sea(l + 7sea) 

Osea—^OO ^TT . , 0 

t=u^d,s 

is given in terms of the normalization of the structure function Fi(a;B, of deeply inelastic 
scattering, see Eq. (|146|) , 


Fi{xb,Q^)^Nf,{Q^)x 


-asea(A2=0,Q2) 

B 


®Since we assume (3 > —1, we can safely set N~ = 0. Note that the pion-pole contribution gives only a subleading 
term in 1/^, however, its numerical value can be rather large for small — A^. 
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Of course, this relation is spoiled by both radiative corrections to the Wilson-coefficients and 
evolution effects. 

In the coefficient of the squared DVCS amplitude the pion-pole contribution induces 

an XB-independent term. As we see from Eq. (|168|), it can be neglected with respect to the 


7i contribution, if the condition ^—/S?/AM‘^ -C 2Ab^/ sin is fulfilled. Indeed, it is 

obeyed in the allowed region of the parameter space for Fgea and Kgea, provided x-q < 0.01 and 
—< 0.5 GeV^. It is interesting that the spin-flip contributions 8 could already be significant 
at low momentum transfer, if /3 > a. Note that further isolated mesonic-like states contribute to 
the small-a:B behavior of For instance, if there would be such a term in 8 with p = 1, we 

will have an additional contribution proportional to -|- ^ 2 ) The effect is even more 

significant in the case of 8, where it would overwhelm the pion-pole contribution and induce an 
interference term with Ti that gives a dominant effect. 

Now let us discuss the extraction of free parameters of our model making use of experimental 
results from HERA. To do this, we narrow down the number of unknowns by assuming the absence 
of isolated mesonic states and consider the small momentum transfer kinematics. Thus, 8 can 
be neglected in the Fourier coefficients (|168|) and ( p.69|) . Then the absolute value of H can be 
extracted from the unpolarized cross section as it has been measured by HI collaboration P] . The 
parameter a can be obtained from the XB-dependence of the measured DVCS cross section. More 
directly it can be gotten from a measurement of the unpolarized azimuthal asymmetries, which 
are proportional to the real part of Ti: 


cos 


TT 

a— I = 


(4,unp 


32K^{2 - 2y + y^)q 


(X 


0,unp 


PpAUnp 


(171) 


Note that for small values of the lepton energy loss y the squared DVCS term is enhanced by 
compared to the squared BH term. Thus, a cleaner separation of the interference and squared 
DVCS term can be done by means of the charge asymmetry. Obviously, the consistency of such 
measurements tests the dominance of both Ti and its sea-quark component, since the measured 
value a should be consistent with deeply inelastic scattering data. However, we cannot exclude 
that a possible A^-dependence of a alters this relation. 

If the polarized lepton beam will be available, one would probe the imaginary part of 
and have so the possibility to measure a directly as a ratio of the real to imaginary part of the 
amplitudes. 


cot Oi— = 


y{2- 2 /)CoA "(“‘5 

(2 - 2|/ + ^ (2 - 27/ -f 7/2 )CoA:-p 


vr\ _ y{2 - y)c{^ 


unp 


(172) 
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Knowing the parameter a, one has an opportunity to access the normalization in any of these 
experiments and, thus, to study the effects of skewedness, which is parametrized by fegea in 
our model. Unfortunately, also the slope i^sea of the form factor F®®®' is unknown, so that the 
skewedness effect cannot be measured in integrated asymmetries. However, due to the extrapola¬ 
tion limAs^o so far it is possible. 

Perhaps, the most interesting quantity is the CFF S since it contains information about the 
orbital angular momentum fraction carried by partons. In our simplistic model the relevant 
parameter is the anomalous magnetic moment Ksea. of sea quarks. Unfortunately, we cannot say 
if S is accessible from unpolarized proton beam experiments at HERA. (A discussion about the 
potential of a polarized proton beam can be found in Ref. m-) From the formulae for the twist- 
two Fourier coefficients (|168|) and ( 169 ) one expects that S can be accessed at larger momentum 
transfer, where, however, the statistics is low. The absence of isolated mesonic states and negligible 
contribution of the spin-flip CFF S allows us to derive a relation between the Fourier coefficients, 

.2 , \ 2 


„DVCS 

■'0,unp 


2-2y + y^ 


2 / 2(2 - yf (4,unp) +i2-2y + y'^f (4,unp)^ 
32F^K^y^{2 - yy{2 - 2y + y^ f 


(173) 


which can serve as a test of our considerations. If the spin-flip contribution S is larger than 
the pion-pole term and compatible with Ti, the deviation from this relation is a measure of Ns. 
Without any further assumptions or theoretical predictions, we cannot directly invert the three 
Fourier coefficients cf, sf, and for an unpolarized target to get the four parameters a, 

N'n, P, and N^. Fortunately, if the quality of data is good enough, one can measure the a,/?- 
parameters by fitting the slope of the xe-dependence. Then again, one has an over-constrained 
set of observables, which can be used to exclude contributions from isolated mesonic-like states. 

Let us now demonstrate how the model parameters can be adjusted by using the HI data 
Pi for unpolarized Compton scattering with the positron beam off the proton, where the hard 
momentum transfer Q and the invariant energy W, 


varied in the range 


W = WM2+ Q‘ 


1 - 
Xb 


(174) 


2 GeV’ < < 20 GeV^, 30 GeV < IP < 120 GeV. 


In this kinematical window the lepton energy loss y is not large, i.e., 0.01 < y < 0.16. The cross 
section (^) is integrated over the momentum transfer squared with the restriction |A2| < 1 GeV^, 
Q2 or W, as specified, and the azimuthal angle is integrated out entirely, 0 < 0 < 27i. For 
simplicity we choose a to be independent on A2, and neglect the evolution and radiative corrections 
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Figure 6: Differential cross sections (a) and ^ (b) for e~^p e'^pj measured by HI 
collaboration [^]. The data are compared with the LO predictions for different parametrizations: 
MRS A’ with bsea oo and Rgea = 9 GeV“^ (solid line) as well as Rgea = 5 GeV“^ (dash-dotted); 
MRS G with Bsea = 9 GeV“^ and bsea oo (dashed) as well as 6sea = 1 (dotted). 


as well as we set Kgea = 0. The integration over the azimuthal angle projects on the constant 
terms in the squared DVGS and interference term. In the hrst case this procedure is exact, 
while for the latter, due to the lepton propagators involved, we only approximately pick out 
the power-suppressed Fourier coefficient cj. Fortunately, for the HI kinematics and within the 
WW-approximation adopted in the estimate, a cancellation occurs between contributions in the 
interference term, and the latter gives a negligible effect (on a percent level) in the total differential 
cross section. The ratio of BH and DVGS cross sections is proportional to (1 — y)\A^\/y‘^Q‘^ and, 
therefore, the BH contribution is much smaller than the DVGS one. In Fig. (a) and (b) we 
show the measured cross sections and respectively, versus model predictions for different 
values of adjustable parameters in our ansatz. Here we take all GFFs into account, however, as 
was explained above, only the sea-quark component of Ti is dominant. 

The measured normalization of the cross section see Fig. (b), can be achieved by tuning 
to a larger value of the slope Bsea, a smaller skewedness effect, or a smaller value of CTsea- The 
best agreement we obtain with the MRS G set of forward parton distributions with Osea = 1.067 
, Bsea = 9 GeV“^, and bsea —^ oc (dashed line in Fig. (|)). Taking the MRS A’ value Osea = 
1.17 we slightly overshoot the data (solid line). The (integrated) DVGS cross section is roughly 
proportional to l/Rgea, thus, for a lower value Rgea = 5 GeV“^ we increase the theoretical estimate 
by almost the factor of 2 for small W (or small y) (dash-dotted). The skewedness effect for the MRS 
G parametrization can be read off from the normalization (|170|) and provides an enhancement of 
the DVGS cross section by the factor of (1.588)^ ~ 2.5 (dotted line). All of these parametrizations 
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Figure 7: Twist-two azimuthal charge asymmetries for e^p —>■ e^pj as functions of for 
VF = 75 GeV and = —0.1 GeV^ (a) as well as functions of W for = 4.5 GeV^ and 
= —0.1 GeV^ (b) for the model of GPDs based on the MRS A’ parametrization with bgea. oo 
and Bsea = 9 GeV“^. For = 2 (/isea = —3) the asymmetries GoA^j'^ and GoA^j'’^ are plotted 
as solid (dash-dotted) and dashed (dotted) lines, respectively. 


are consistent with the measured differential cross section Fig. (^ (a). The evolution flow 
with Q^, which we presently neglect, will give a moderate logarithmic change of H with increasing 
Q^. Note that according to Ref. H the data can also be £t in the aligned jet model and by 


adopting the reggeon and soft pomeron exchange mechanism ^0 


Since the DVGS amplitudes enter the DVGS cross section in square, the size of NLO corrections 
there is twice of what we have found for the GFF Ti in section |6.2.3| . They are of order —40% even 
in the DVGS scheme. This requires a NLO analysis in order to extract the NLO parameters of the 
GPD models. However, since our ultimate goal is to compare twist-two and -three approximation, 
with the latter available only in LO of perturbation theory, a consistent approach is to use the LO 
parameters we have extracted. We will stick in the following to the MRS A’ set with Aea —^ cxo 
and Bsea = 9 GeV“^. Note that the dependence of predictions on the Kgea parameter is weak and 
was previously neglected. However, in the following we can use it also as a free fitting parameter. 

As in the case of cross sections, the size of asymmetries can easily be estimated from the 
analytical formulae we have presented. The advantage is obvious: (i) one gets an understanding of 
the size of asymmetries depending on kinematical changes and modihcation of model parameters, 
(ii) one can easily scan the parameter space quantitatively. 

The magnitudes of the twist-two azimuthal charge asymmetries GoA^j*^ and GoA“|)^ are dis¬ 
played in Fig. ^ for A^ = —0.1 GeV^. The absolute values of the asymmetries arising from the 
cos(0)- and sin(0)-dependence in the interference term can reach almost 15% and 10%, respec- 
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Figure 8: Twist-three azimuthal charge asymmetries for the same settings as in Fig. (a) and (b), 
respectively, are plotted for GeA^j*^ (solid), GoA^^q*^ (dash-dotted), (dashed), and GoA™ 2 ’) 

(dotted). 


tively. It is interesting that the asymmetries are reduced by about 30% if we vary the parameter 
Ksea from 2 to —3. This dependence arises mainly from the interference term, where the S- 
contribution is enhanced by the normalization of F 2 , i.e., by the anomalous magnetic moment of 
the proton [cf. Eqs. (|169|) ]. Gonsequently, even with a low mean value of the momentum transfer, 
e.g., —0.1 GeV^, the T-contribution can be signihcant for charge and lepton single-spin 

asymmetries. Note that for Ksea = —3 in our ansatz, T contributes already for —0.6 GeV^ 

with the same magnitude as —7^. Asymmetries, averaged over the momentum transfer with the 
restriction |A^| < 1 GeV^, have almost the same size as the ones displayed in Fig. |^, e.g., 

GoA“P(fF = 75 GeV, = 4.5 GeV^) ^ 8%, 

GoA“P(hF = 75 GeV, = 4.5 GeV^) ^ -2%, 


however, with a weak Kgea-dependence. 

Now we want to answer the second question we raised at the beginning of this section about 
the size of twist-three effects. In the WW-approximation for our model only the sea-quark part 
of 7i is dominant. Moreover, from Eqs. (|84[) and (|87D we see that 7-f-*“(^, A^) drops out for small 
x-Q- Thus, in the limit 6sea —cxo we discover from the WW-approximation of Eq. (Q) with the 
help of Eqs. ( |157]) and ( |161| ), 


n 


eff-WW 


(e,A2) = 2-a 




1 + dg 




Qfc 


WK,A5, 


(175) 
= 1.17 the 


2 J " V 2 

a known relation, which was derived in the archive version of Ref. [Q. For CKg^ 
ratio is ~ 0.65, while for hsea = 1 we hnd ~ 0.48. Apparently, the size of the twist- 

three effects is entirely given by the kinematical factors in front of the ‘universal’ C functions in 
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model 

^val 

^sea 

Rsea [GeV-2] 

^sea 

A 

1 

OO 

9 

0 

B 

OO 

OO 

9 

-3 

G 

1 

1 

5 

0 


Table 3: Parameter sets for models of H- and T^-type GPDs. 


the definition of Fonrier coefficients (|53|) , and (|55|) . Thus, we realize that the twist-three 
harmonics are relatively suppressed by the factor \J—{1 — y)^‘±/Q^- For instance, the suppression 
is already about the factor of 4 for = —0.1 GeV^ and = 2 GeV^. In Fig. ^ we demonstrate 
that this counting is not justihed however for the asymmetry GoA“^^^ as one can easily realize from 
Eq. (^). For = —0.1 GeV^ we hnd a contribution (dash-dotted line) that is of order of a few 
percent. Fortunately, it is expressed in terms of twist-two GFFs and so this twist-three azimuthal 
asymmetry can be employed to pin down twist-two GFFs. Within our model we expect that the 
addendum (l7^) does not contribute in Eq. (^) and, thus, we have the ratio 


CoAqT.'S ^ {2-y? -Aq ( AL„\ 2(1-!/) 1 

CoA"”!; 2(2-2y + y^)Q^K\{ ) (2 - ' 


The largest size among the twist-three azimuthal asymmetries has GeA^Ij’j (solid line), which 
arises from the cos(0) term of the squared DVGS term. Its twist-two counterpart GeA^^'IJj ap¬ 
proaches the value 1 for small W (small y), which reflects the fact that the DVGS process domi¬ 
nates over BH one. If we average over the momentum transfer squared, we still And a few percent 
effect for the twist-three azimuthal charge asymmetries that is comparable with the twist-two 
ones. 


GeA"“P(IE = 75 GeV, = 4.5 GeV^) ^ 7%, 

GoA“P(iy = 75 GeV, = 4.5 GeV^) ^ 2%, 

for -AT^ < |A2| < 1 GeVl 

7.2 Asymmetries in unpolarized fixed target experiments 

After fixing the parameters in the sea-quark sector from the HI data, we could hope to give better 
numerical estimates also for unpolarized fixed target experiments. To have a comprehensive 
understanding of this issue, we define three parameter sets for H and £ based on the MRS A’ 
parametrization. They are summarized in Table The slope parameter Bd = S/m^ of the 
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D-term will be fixed later. In Ti we use the parameters 


&vai = &sea = 1, Ba = 3/m^ Ri 3.7 GeV 


-2 


Obviously, in the models A and B, which are compatible with the HI measurements, the sea- 
quark contributions are relatively small. For the model C this is not the case and our predictions 
overshoot the HI data. If our oversimplified model has some realistic features, the model C 
will hardly describe the beam-spin asymmetries measured at HERMES and Jefferson Lab. The 
philosophy would be then again to employ experiments to fix the model parameters left in the 
valence-quark sector or to discriminate between the model A and B, or any other. 


In the following section [7.2.1| we compare these models, used for estimates of the DVCS am¬ 
plitudes calculated in the twist-three approximation, with the experimental measurements and 
discuss theoretical uncertainties in the extraction of model-dependent parameters. Then we have 
in section [7.2.2| a closer look to the twist-three effects. Especially, we show that asymmetries, 
usually discussed in the literature, are not quite suitable to access GPDs. In the third section 
[f.3| we estimate the size and parameter dependence of nucleon-spin asymmetries for the GLAS 
kinematics. 


7.2.1 Beam-spin asymmetry measured with HERMES and CL AS 

The azimuthal angular dependence of the beam-spin asymmetry, i.e., 

da\(l)) - da^^cj)) 

has been recently measured in two different fixed target experiments, namely, by HERMES and 
GLAS [Q collaborations with the positron and electron beam of 27.6 GeV and 4.25 GeV, respec¬ 
tively, scattered on hydrogen targets. At HERMES the average values of kinematical variables are 
(Q 2 ) ^ 2.6 GeV^, {x-q) = 0.11 and (—A^) = 0.27 GeV^ |^. The results at GLAS are integrated 
over the regions 1 GeV^ < < 1.75 GeV^, 0.13 < xb < 0.35, and 0.1 GeV^ < —< 0.3 GeV^ 

with the condition that IV > 2 GeV 

In both experiments we have a large average value of (y), i.e., ~ 0.5 and ~ 0.85 for the 
HERMES and GLAS measurements, respectively. Thus, we expect that the BH part of the 
total cross section dominates over the DVGS one [like (1 — y)A‘^/y'^Q^ as indicated by Eqs. 


^)], which is indeed the case. So the denominator, da^ + da^, of the beam-spin asymmetry 
is prevailed by the BH contribution with Cq ^^p-coefficient being the dominant one and 
higher harmonics suppressed by one or two powers of K. The BH part provides a much stronger 
contribution for the GLAS than the for HERMES experiment. Thus, we realize that the 0- 
dependence of the BH propagators in Alu( 0) almost cancels. From this qualitative discussion 
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Figure 9: The beam-spin asymmetry as function of the azimuthal angle measured in e^p e'^p'j 
with E = 27.6 GeV from HERMES (a) and e~p e~p'j with E = 4.25 GeV from GLAS (b) 
collaborations. The dotted lines show the sin(0) function with the amplitude 0.23 and 0.202 for 
HERMES and GLAS, respectively. The other curves are predicted by the model A (solid) and 
G (dashed) in the WW approximation as well as the model B beyond this approximation (dash- 
dotted), as specihed in the text, for = 2.6 GeV^, xb = 0.11, and — = 0.27 GeV^ in (a) as 
well as for = 1.31 GeV^, x-q = 0.19, and —= 0.15 GeV^ in (b). 

we expect that in the kinematics, we are considering, the azimuthal dependence is approximately 
given by the sin(0) function: 

Alu(0) ~±—%^sin(0) oc Am I E^H + (Fi + E^)^ - | sin(0), (177) 

y Co,unp ( ^ J 

where higher harmonics are suppressed by A/Q. Moreover, Alu(0) in this approximation is linear 
in GFFs, and in the kinematics, we are discussing, the dominant contribution arises from AmT^ 
within our model. In LO the imaginary part is directly given by the GPDs on the diagonal 
X = ±.^. Thus, uncertainties due to our poor knowledge of GPDs are minimized for the beam-spin 
asymmetry. Moreover, after hxing the parameters for TYsea from the HI data, we can confront 
our predictions with the data. Obviously, the D-term can only slightly affect the beam-spin 
asymmetry. Nevertheless, we include it and set the sea- and D-slope parameters equal Bjy = i?sea- 
Now let us compare our expectations at the twist-three level for LO GFFs with the results 
from the two experiments shown in Fig. p. 

The model predictions for the GLAS experiment, plotted in Fig. (b), possess a much smaller 
deviation from the sin(0)-like shape than the HERMES one shown on the panel (a). For this 
kinematics a stronger contamination by higher harmonics, due to the higher average (?/), can be 
induced by the squared DVGS term in da^ -|- da^ as well as twist-three effects. The model A is 
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compatible with the HERMES data (solid curve). For the integrated asymmetry 




7 


the experimental result from HERMES 


Hlu = —0.23 ± 0.04(stat) ± 0.03(syst), 


is in agreement with the model-A prediction 


Hlu = -0.27. 


( 178 ) 


If we take a smaller slope parameter Rgga or a smaller 6gea parameter, the absolute value of the 
asymmetry will increase. For instance, the model C leads to 


Hlu = -0.37 


with an obvious deviation from a simple sin((/))-shape (dashed), which is mainly caused by purely 
kinematical effects. This parameter set yields a prediction that is close to the result of Ref. . 
On the other hand we can also decrease the absolute value of Hlu if we take a larger value of 6vai, 
a smaller Ksea., or if we go beyond the WW-approximation. To demonstrate possible effects due 
to the antiquark-gluon-quark correlations, we assume a rather large contribution proportional to 
the twist-two one, however, with a unique phase difference]^ 

jrqGq ^ 1 |jr| exp {i arg(J^) + , for = {h, S, h] . (179) 

Here we also included an extra factor 1/,^ since the convolution given in Eq. (|8^) looks rather 
singular. Now the model B with = —n/S results in a smaller asymmetry, 

Hlu = —0.16, 


which is also compatible with the data (dash-dotted curve). In this case, we observe the appearance 
of sin(20) and sin(30) harmonics, which arise from the interplay of twist-three effects in the 
interference term and purely kinematical effects. Note that the amplitude of the sin((/))-harmonic in 
Hlu( 0) is insensitive to the twist-three sector]^. However, for the azimuthal beam-spin asymmetry 
(|178|) this does not necessarily hold true. 


have to emphasize that this is only a toy model in which possible constraints due to sum rules and reduction 
to the forward kinematics for antiquark-gluon-quark GPDs were not implemented. 

^^Such dependence arises mainly from twist-three effects in the denominator on the right-hand side of Eq. (176), 
which is, however, dominated by the BH cross section. 
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Now we come to the measurements of CLAS collaboration, which can be £t by 


^Lu(0^y) = asin(0^y) + /3sin(20^^), 


withQa 


a = 0.202 ± 0.028(stat) ± 0.013(syst), 
f3 = -0.024 ± 0.021(stat) ± 0.009(syst). 


The /3-parameter is obviously compatible with zero. In the case of the model C (dashed curve) 
and now also for the model A (solid curve) one fails to describe the data since one gets 


a = 0.4, /3 = 0.028, 


and 


a = 0.28, /3 = 0.014, 

respectively. Taking now the model B and going beyond the WW-approximation by inclusion of 
antiquark-gluon-quark effects with the phase difference = vr, we find a result compatible with 
the data (dash-dotted curve), namely, 

a = 0.24, /3 = -0.03. 

As we mentioned before, the value of a is almost insensitive to the antiquark-gluon-quark contri¬ 
bution. Indeed, in the WW-approximation we have for the model B the same value of a but quite 
different /3, i.e.. 


a = 0.24, /3 = 0.014, 

which is slightly out of the experimental range. Our model-dependent analysis shows that the 
CLAS measurement of the beam-spin asymmetry is sensitive to antiquark-gluon-quark contribu¬ 
tions. However, due to the large experimental uncertainties, this statement might be premature. 
On the theoretical side, the legitimacy of this assertion depends on the answers to the following 
two crucial questions: (i) Can we be sure that our equations are valid at low as of order of 1 
GeV^? (ii) If yes, is the sign of f3 model-independent in the WW-approximation? 

To address the first question, we give now a more qualitative discussion about our present 
understanding of corrections, which are expected to modify our predictions. We provide here 
not mix a and /3 in this parametrization with y and ^ exponents of parton densities and GPDs, respectively, 
discussed in the preceding sections. 
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Figure 10: Estimates of the beam spin asymmetry 74 lu( 0) for the E-91-023 experiment at 
Jefferson Lab with the E = 6 GeV electron beam scattered on a proton target. Predictions for 
the model A (solid) in the WW-approximation and B (dash-dotted) in the same approximation 
as well as with the antiquark-gluon-quark contribution with = tt (dashed) and = 0 
(dotted) are plotted in (a) for xb = 0.3, = —0.25, and Q = 2.5 GeV^. Here the D-term has 

been taken into account with Bd = Hsea- In (b) we show the model B estimates with = tt for 
Xb = 0.35, = 3 GeV^ (solid), xb = 0.3, = 2.5 GeV^ (dash-dotted), xb = 0.25, = 2 GeV^ 
(dashed), and xb = 0.19, = 2 GeV^ (dotted), where in all cases = —0.25 GeV^. 


some hints rather than answers. In both experiments the perturbative NLO effects in the DVGS 
scheme may be of order of —20% or so (cf. Fig. ^). Furthermore, we expect that the evolution 
effects can alter predictions on a few percent level at most. The power corrections involved in 
the analysis are of two kinds —A^/and M^/ Q^. While the hrst correction in our model is 
naively expected to be of ±10% or so in both experiments, the latter is certainly larger for the 
GLAS settings and, therefore, their quantitative estimate is desirable. To this end, one can adopt 
the partial result of Ref. [|^. In contrast to the forward case, where each power of the target 
mass is accompanied by a suppression factor x|, the situation in the DVGS kinematics is 
not so obvious and requires detailed numerical studies. So far it is likely that besides radiative 
corrections also power-suppressed contributions are important for the interpretation of the GLAS 
experiment. This is especially crucial for the interpretation of the sin(20) harmonic. 

Keeping these reservations in mind, we hnally estimate the beam-spin asymmetry Alu( 0) that 
we expect from future measurements at Jefferson Lab with the electron beam energy E = 6 GeV 
in E-91-023 experiment with GLAS. Our expectations are shown in Fig. (a) and (b) for the 
models A and B in different approximations as well as for different average values of the kinematical 
variables. They are chosen so, that 0.7 < y < 0.76 and, thus, the DVGS amplitude is kinematically 
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suppressed. We have a similar situation as discussed above, however, for larger Q^. This may 
offer a possibility to investigate a deviation from the predicted asymmetries, expected in the 
WW-approximation, in more detail. To clarify the issue whether this deviation is caused by the 
antiquark-gluon-quark or the mass corrections, it is necessary to £x xb and vary Q^. For instance, 
our expectation for the average value of x-q = 0.19 (as in Fig. (b)), are shown as the dotted 
line in Fig. 0 (b) for = 2 GeV^. For this kinematics, we have only a tiny deviation from the 
sin(0)-shape in the WW-approximation, i.e., 

^lu( 07) = 0.136 sin(0!^) -F 0.005 sin(20(y). 

A shape closer to a sin(0)-like form in the E-91-023 measurement would indicate that the source 
of higher harmonics in the CLAS experiment at lower beam energy is caused by mass corrections, 
since they are now relatively suppressed by the factor of ~ 0.65. Obviously, a dehnite statement 
requires higher statistics in order to decrease the error on jS. 

7.2.2 A closer inspection of twist-three effects 

The essential feature of the twist-three contributions is the appearance of new coefficients including 
a 0-independent harmonic in the interference term. It can be seen from the equations given in sec¬ 
tions and that for longitudinally polarized target higher Fourier coefficients in the squared 
DVCS as well as interference term are relatively suppressed by the factor K oc ^(1 — |/)A ^/Q 
with an additional multiplicative ^/-dependence. Thus, higher harmonics vanish faster at the kine- 
matical boundaries —>■ 0 and // —>■ 1 than the twist-two ones. The nonperturbative dependence 

of the twist-three CFFs is given by the universal C coefficients. In the WW-approximation these 
CFFs have a magnitude similar to the twist-two ones (see Fig. |^). Although for the longitudinally 
polarized target the constant interference terms cj, defined in Eqs. (|5^) and (^71) , are suppressed 
by 1/Q, they do not vanish at the kinematical boundaries. Therefore, we expect the existence of 
phase space regions where the interference term is strongly contaminated by this purely kinemat¬ 
ical twist-three effect. Note that this is not the case for cj^p, see Eq. (^). 

Let us study the twist-three effects for the longitudinally polarized fixed target experiment in 
more detail. For = —0.25 GeV^ and = A GeV^, we expect a suppression of the higher 
harmonics, arising beyond the twist-two approximation, roughly by the factor of 4. Note that 
A^/Q^ ~ 0.06 and the target mass effects, i.e., could be suppressed by a similar factor 

since they are combined with an extra multiplicative ^-dependence. So in the future precise DVGS 
measurements there will hopefully be a kinematical window left, where twist-three effects are 
accessible and 1/Q^-power corrections are quantitatively not important. We set for the following 
numerical estimates Ebeam = 27.6 GeV, and take the model A with the pion pole and D-term, 
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having equal slopes = -Bsea- For the squared DVCS term we hud, indeed, the anticipated 
suppression in the WW-approximation. For the longitudinally polarized target and polarized 
lepton beam, the Fourier series for xb = 0.1 (?/ ~ 0.77) reads 


dx-Bd\A‘^\dQ^d(j) 


0.0023 [1 + 0.19 cos(0) - 0.09A sin 


( 180 ) 


f nb 


+0.0004A [A + 0.16A cos(0) — O.O4sin(0)] j 

Let us comment on these numbers. The normalization of the cross section strongly depends on 
the parameter set we have chosen. For instance, for the model C we find an enhancement in the 
normalization almost by the factor of 3. However, it depends only weakly on the pion-pole and 
D-term contributions. The suppression of terms involving target polarization arises rather from 
the destructive interference of Ti and Ti and should be considered as a reflection of the model 
taken for the estimate. For the cos(0) harmonics we have the expected damping, while the strong 
suppression of the sin((/)) terms arises from the adopted WW-approximation. The main ingredient 
of the latter leading to this feature is the fact that the twist-three CFFs have a phase similar 
to the twist-two ones. Thus, the imaginary part of their interference is roughly zero. A strong 
deviation from such a small value would indicate an antiquark-gluon-quark contribution having 
quite a different phase structure as compared to the twist-two CFFs. 

The interference term for the positron scattering reads 

da^ 1 


X 


1 -|- 0.93 cos(0) -|- 0.08 cos(20) 

0.007 [0.56 + cos(0) + 0.17 cos(20)] - 0.02A [sin(0) + 0.06 sin(20)] 
+0.002AA [0.52 + cos((/)) + O.O3cos(20)] - 0.006A [sin(0) + O.O9sin(20)] 


( 181 ) 


nb 


GeV^ ’ 

with the same kinematical settings as above. Again, the overall normalization strongly depends on 
the GPD models. The suppression of higher harmonics is of the magnitude we have expected or 
even larger, by the same token as before. The constant term is rather large. As already mentioned, 
this is a purely kinematical twist-three term, i.e., it is completely determined by the twist-two 
GFFs. So we see that such a term can alter the cross section in a much stronger manner than the 
others. However, it will not spoil the separation of the dynamical twist-two and -three sectors. 
Thus, it is helpful for both experimental consistency checks and access to the real part of GFFs. 
Let us also note that the 0-dependence arising from the BH-propagators is rather strong for the 
kinematics, we have chosen. 


It is instructive to compare our results in the WW-approximation with the one given in [^ 
The sets of GPDs used there presently differ by the choice of the forward parton densities and 
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da/{dQ‘^dA^dxBd(j)) [nb/GeV^] da/{dQ‘^dA^dxBd4>) [nb/GeV^] 





Figure 11: The unpolarized e~p —^ e~pj cross section (a), the charge asymmetry (b) for = 
2.5 GeV^, = —0.25 GeV^ and x-q = 0.3 as well as the polarized electron-beam cross section (c), 
and the single-spin asymmetry (d) for = 2.5 GeV^, = —0.25 GeV^ and xb = 0.15 are shown 
for the BH term (dotted), the twist-two without the D-term (dashed), twist-two (dash-dotted) 
and -three (solid) with the D-term for electron beam within the model G. 


partonic form factors. To be close to the settings of Ref. [^, we take the model G with Bb, = 
Bsea. = 5 GeV“^. Note that our Hgea increases less rapidly at small ^ than in the parametrization 
from Ref. [^, and we choose a stronger fall-off with A^ of the D-term. In Fig. |^, we show 
the same predictions for the (/)-dependence of the cross section as given in Figs. 2 and 4 of Ref. 


ISj . Our estimates for the unpolarized cross section (a) and the polarized beam (A = 1) cross 
section (b) for electron scattering on the unpolarized proton target with = 2.5 GeV^ and 
A^ = —0.25 GeV^ look quite similar to IM. We can clearly observe large differences of the cross 


sections and charge asymmetry, i.e., {da~^{(j)) — da~{(p)}/{da~^{(j)) + da~{(j))}, between twist-two 
(dash-dotted curve) and -three (solid curve) contributions. As we have shown previously in Eq. 
(|181|) , the main effect comes here from the interplay of the constant interference term with the 
^-dependence of the BH propagators. 


62 





















The most significant differences of both resnlts are prominent in the charge azimnthal asym¬ 
metry, displayed in Fig. |TI| (b). Since in our parametrization the positive sea-quark contribution is 
larger, we hnd without the D-term a large negative value of this asymmetry at = 0 in the twist- 
two approximation (dashed curve). Adding the D-term contribution we get the asymmetry, which 
is compatible with zero in the same approximation (dash-dotted curve), while the parametrization 
of [jT^ yields a large positive value at 0!) = 0. So we conclude that due to many parameters involved 
in the modeling of GPDs there is no unequivocal experimental signature of the D-term. Moreover, 
at the twist-three level without the D-term we have a similar magnitude of the asymmetry as 
displayed by the solid curve in Fig. (b), which includes the D-term. This tells us that there is a 
large cancellation between the twist-two and -three effects induced by the D-term. Consequently, 
to access the twist-two GPDs one has to separate the twist-two and -three sectors. Note that the 
cancellation between sea- and valence-quarks as well as the D-term contributions induces a rather 
strong dependence on the parametrization for the unpolarized beam-charge asymmetry, and, thus, 
may also induce a stronger contamination by radiative corrections as was observed for separate 
contributions. 

The beam-spin asymmetry in Fig. 0 (d) is only slightly modihed by the twist-three contri¬ 
butions in the WW-approximation. These corrections enter mainly in the denominator of this 
asymmetry, while the contamination of the numerator by new harmonics remains small as we 
have discussed above in Eqs. (|18CI|) and (|181|). Again, for this asymmetry we expect only moder¬ 


ate changes due to radiative corrections within our ansatz of order of 20% or so. 


7.2.3 Estimates and parameter dependence of charge asymmetries 

As we have demonstrated in the previous two paragraphs for the beam-spin and charge asymme¬ 
tries, an important phenomenological issue is the separation of the twist-two and -three sectors. 
This problem is completely resolved by the introduction of charge asymmetries via Eqs. (|105| - |109|) 
and partially for the single-spin asymmetry defined by Eq. O- We emphasize again that the 
charge asymmetries allow also to access the imaginary part of the interference term in a clean 
manner, which is hardly possible for the beam-spin asymmetries. Differences encountered due to 
various definitions of the charge and beam-spin asymmetries in the twist-two and -three approx¬ 
imations are shown in Fig. 0 - The asymmetries Ac (a), Asl and SSAi (b) are normalized to 
the unpolarized cross section that does not or only weakly does depend on the dynamical twist- 
three contributions. However, as was pointed out above, the constant term in the interference 
may generate an important effect. From the differences of the twist-three (dotted) and twist-two 
(dash-dotted) approximation one observes in the panel (b) that the twist-three contamination is 
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Figure 12: The azimuthal charge asymmetry Ac (a) and the beam-spin asymmetry Asl (b) 
are plotted in twist-two (dash-dotted) and twist-three (dotted) for the HERMES kinematics with 
= 2.5 GeV^ and = —0.25. The solid curves show the predictions of the improved charge 
asymmetries CoA“^”|^ and CoA^I)*^. In (b) the short-dashed and dashed curves display SSAi in twist- 
two and -three approximation, respectively. All curves are given for model A with Rd = 9 GeV^ 
and includes an antiquark-gluon-quark correlation with = —n/3. 


quite essential for smaller values of xb, i.e., larger value of y, where the BH propagator Vi{(f)) 
induces a rather strong 0-dependence. In the case of the beam-spin asymmetry an improvement 
is necessary if considerable antiquark-gluon-quark correlations are present. As it can be seen by 
comparison of the twist-two (short-dashed) and -three (dashed) approximation for SSAi, a clean 
separation is only achieved when the DVGS amplitude is suppressed with respect to the BH one, 
i.e., \J—{1 — y)/SP‘/yQ -C 1. Note also that the twist-two approximation of SSAi is close to the 
charge asymmetry GoA“|)^. The deviation of GoA^I)*^ from SSAi arises only from the interference 
term in the denominator. 

Unfortunately, the improved dehnitions reduce the size of asymmetries. This disadvantage is 
unavoidable provided one is willing to extract in a clean way information on the twist-two and 
-three GPDs. 

Now let us discuss in more detail the magnitude of the improved charge azimuthal asymmetries, 
dehned in Eqs. (|lU5[ - p.l2| ) for the charge-odd and -even parts of the cross section, in dependence on 
the twist-two parametrizations and antiquark-gluon-quark correlations. For the twist-two sector, 
we present our estimates in Fig. with the models A and G. As was mentioned above in section 


7.2.2| , in contrast to the sin(0) harmonic the size and the sign of the cos(0) term crucially depends 
on the parametrization of the valence- and sea-quark densities as well as the D-term contribution. 
According to our model, the sea quarks give a positive real part of H, while the valence quarks 
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Figure 13: Azimuthal charge asymmetries of cos-harmonics CoA™^^ (a), CeA“^Q^ (b), CoA“|^^ (c), 
and CoA™^ (d) are plotted for i^beam = 27.6 GeV, = —0.25, and Q'^ = 2.5 GeV^. The model 
A (G) with the D-term, Bb, = 5 GeV“^, and neglected one are displayed as dash-dotted (dotted) 
and solid (dashed) curves, respectively. 


and D-term provide a negative effect. The GA°^]^ plotted in panel (a), is now positive if it is 
dominated by the sea quarks, like in the model G without the D-term (dashed), contrary to all 
other cases. There is also a possibility that those contributions may compensate each other and 
this will lead to the almost vanishing asymmetry. Since the model G (dashed) is ruled out by the 
beam-spin measurement, we rather expect a sign alternating asymmetry [the model A without the 
D-term (solid)] or a negative one [the model A with Bb = 5 GeV“^ (dash-dotted)]. Its magnitude 
can be of order of ±10%, however, it can also be compatible with zero. In panel (b) we plot the 
leading twist-two asymmetry GeA]][)Jj, which shows that the BH amplitude dominates for small 
Xb- Also for HERMES there is a kinematical region in which the DVGS cross section overwhelms 
the BH one. This happens for the model A (B) at xb ~ 0.2 [xb ~ 0.17). The DVGS cross section 
only weakly depends on the D- and pion-pole term. 

In Fig. (c) we show the twist-three asymmetry GoAc(o),unp, which can be rather sizeable, 
especially, at the kinematical boundary 0. The complexity of its shape partly emerges 
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Figure 14: Twist-three azimuthal charge asymmetries of cos-harmonics CoA“|^^ (a) and CeA^j*^ 
(b) are plotted for the same kinematics as in Fig. |^. The model A with (without) the D-term 
in the WW-approximation and with the antiquark-gluon-quark correlation = —vr/S) are 

displayed as dash-dotted (dotted) and solid (dashed) curves, respectively. 


from the kinematical factors in Eq. (^). According to this equation and from the dehnition of 
the charge asymmetries in section b.2.1|, we have the constraint 


CoAr = CoA""P 


TT (1 -XB){2-yf -A^ I 


TT 


2 (2 
(2 - y) 


2(2-2y + 2 / 2 ) 


2y + 2 / 2 ) 

(1-2/)(2-xb) X 


^Ln \ _ ( 2 - 3 : b )( 1 - 2 /) 

A2 ) (1 -a;B)(2-2/)2 


>CoA-P 


(182) 


where the ellipsis stand for terms generated by the addendum AC;f^p, dehned in Eq. ([7^). For 
large y this difference should vanish. Although, ACy^p is suppressed by at least one power of xb, 
we see in Fig. |T^ (d) that it starts to be sizeable for xb > 0.3. This effect arises from the fact that 
for this kinematics in our model \H\ > |7-f|, and Ti can contribute considerably both to CA°)\^pj^p 
and CoA))^“(J^. If one would measure a large deviation of CoA™^ from zero, it will indicate that 
3?e7i is small and does not give the dominant contributions in the even harmonics. 

Now we come to the dynamical twist-three effects, which show up in the charge asymmetries 
CoA))^ 2 *) CeA™]*^ displayed in Fig. |T^ (a) and (b), respectively. In the WW-approximation with 

the model A (dash-dotted and solid curves) we hnd again the expected magnitude of asymmetries: 
the size of the charge-even asymmetry CeA“^”|j is larger than the twist-two charge-odd one CoA^j'j. 
Furthermore, we see that the shape of CoA))^^*) is correlated with the one of CoA™^^. A large CoA))^ 2 *) 
asymmetry can be generated by multi-parton correlation effects shown as dashed and dotted lines, 
respectively. Such correlation can even enhance the asymmetry CeA^^^^. 

Let us add that the parameter dependence in the twist-two sector and twist-three effects for 
the polarized beam can also be discussed in a clear manner in terms of the charge asymmetries 
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Figure 15: Estimates of the target spin asymmetry 74 ul( 0) for the E-91-023 experiment at 


Jefferson Lab [|I3 with the E = 6 GeV electron beam. Same models and kinematics as in Fig. |l0 
(a) and (b). 


CoA^j’’^, CoA™ 2 p and CeA“|j'j. However, since qualitatively they display the same effects as it 
has already been discussed for the beam-spin asymmetry in section |7.2.1| , we will not repeat these 
considerations here. 


7.3 Numerical estimates for longitudinally polarized target 


Finally, we provide quantitative estimates for the cross sections of polarized electron scattering 
off the longitudinally polarized proton target. In these settings there are two further observables, 
which can be used to unravel GPDs from data. 

The hrst one is the target-spin asymmetry with unpolarized lepton beam 


da'^{(j)) + da^{(j)) ’ 


(183) 


where the up (down) arrow stands for A = -|-(—)1 longitudinal polarization. Gompletely analo¬ 
gous to the beam-spin asymmetry, the BH cross section cancels and one is left with the leading- 
and higher-twist contributions from the interference term and power-suppressed effects from the 
squared DVGS amplitude. From the explicit expressions for the Fourier coefficients one can qual¬ 
itatively understand the magnitude of the asymmetry. For instance, comparing the main term 
in the asymmetry C^p, see Eq. with Cunp, previously analyzed in section |7.2.1|, one observes 
that the GFF Td, dominating the latter at moderate and small xb, now enters the amplitude with 
an additional power of xb- It becomes parametrically of the same order as the parity-odd GFF H: 
\H\ ~ XB|7f|. Thus both of them play a distinctive role in building up the nucleon-spin asymmetry 
which is displayed in Fig. |^for the kinematics of the E-91-023 experiment at Jefferson Lab with 
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Figure 16: Estimates of the double spin asymmetry y4LL(0) for the E-91-023 experiment at 


Jefferson Lab [|T3 with a E = 6 GeV electron beam. Predictions for the model A (solid) with 
i?D = Bsea in the WW-approximation and model B in the same approximation without D-term 
(dash-dotted) as well as with D-term and antiquark-gluon-quark contributions with = tt 
( dashed) are plotted in (a) for xb = 0.3, = —0.25, and Q = 2.5 GeV^. The dotted curve shows 

the BH contribution alone. In (b) we show the model A estimates for the same kinematics as in 

Fig. 0(b). 


E = 6 GeV electron beam |W|. On the left panel (a) we give the estimates with the models A 
and B for the kinematical variables chosen as xb = 0.3, = —0.25 GeV^, and = 2.5 GeV^. 

While on the panel (b) we display the asymmetry for several settings of the kinematical variables 
with the model B, which gave the best agreement with the unpolarized proton data. 

The second interesting observable is the double lepton-nucleon spin asymmetry. The combi¬ 
nation of the cross section which picks out only the AA-part, bilinear in the nucleon and lepton 
helicity, is achieved by the following combination of the leptoproduction cross sections 

(Jcr^^(0) — — da^^{(j)) + da^^{(j)) 


All(0) = 


(J(jfo(0) -|- dcrfo(0) -|- da^^{(j)) + (JcrD(0) 


(184) 


Since the double-spin asymmetry is not protected from the contribution of the BH process, one 
expects that it overwhelms the DVGS signal inasmuch as it dominates the cross section for the 


GLAS kinematics. Indeed, this is what one hnds in Fig. ^ (a), where the BH cross section alone 
generates a large asymmetry^, displayed as dotted curve. Obviously, the interference and squared 
DVGS term produce a signihcant signal on its background. Due to rather good knowledge of the 
BH process its subtraction from the data should not present an obstacle of principle. However, 
what is more troublesome is the contamination of the asymmetry by the DVGS cross section and 


^^Note that we set in the asymmetry All the DVGS amplitude in both the numerator and denominator to zero. 

















the fact that three CFFs, i.e., 7i, 7i, and the pion-pole term of 8^ can contribnte on the same 
footing to the interference term. This makes the disentanglement of separate components and, 
therefore, GPDs unmanageable from such an experiment alone. The expectations for the same 
kinematical settings as in Fig. (b) are given in Fig. (b) for the model A. 


8 Summary and conclusions 

Our present study has resulted into a complete analytical structure of the real-photon leptopro- 
duction cross section with power accuracy. Although these results are quite complex, they allow 
for a qualitative discussion of observables in dependence of the kinematics and GPD parameters. 
Thus, in contrast to a purely numerical approach, our analysis tolerates a careful examination of 
experimental data. The most accurate information on GPDs can be deduced from the part of the 
total cross section, which stems from the interference of the BH process with the DVGS amplitude. 
The main tool of experimental exploration of the latter is the use of diverse asymmetries involving 
charge and spin. 

The leading contributions in the light-cone expansion of the DVGS amplitude give access to 
diverse properties of the nucleon constituents due to a non-zero orbital momentum carried by them. 
These include, the parton angular momentum, the helicity-flip polarized glue uncontaminated by 
quarks, the pion cloud of the proton etc., just to name a few. 

Twist expansion of the DVGS amplitude hnds its reflection in the dependence of the cross 
section on the azimuthal angle of the outgoing real photon (or recoiled proton). One of the main 
results of our analysis is the hnding that the Fourier harmonics of the separate charge-even and 
-odd components of the complete cross section, expressed in terms of the lowest-twist GPDs, stay 
uncontaminated even if 1/ Q-power effects are switched on. Thus, power corrections to the Fourier 
coefficients will be suppressed by 1/Q^. This gives a hope for a clean exploration of individual 
contributions, provided the same settings will become available in experiments. This can be 
achieved by a combining use of charge and spin asymmetries together with a precise extraction of 
different Fourier coefficients. These procedures enable a separate measurement of all coefficients 
of the interference and squared DVGS term. Extraction of the latter requires the subtraction of 
the squared BH amplitude from data. 

For a lepton beam with a dehnite charge, high precision data may also result into a separation of 
the interference and squared DVGS terms made feasible by a Fourier analysis. Albeit, in this case 
the leading twist cos(0)/sin(0) dependence in the cross section, stemming from the interference 
term, gets corrected by the twist-three contribution of the squared DVGS amplitude. 

The magnitude of higher-order corrections strongly depends on the models assumed for GPDs. 
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The minimization of the radiative effects is possible within the so-called DVCS scheme. It is 
dehned by the condition of absorption of gluonic higher-order corrections into the factorization 
scale, dividing short and long distances. This procedure results into a moderate modihcation of LO 
predictions, advocating the neglect of NLO terms for a crude estimate. Once a more accurate data 
will become available, a NLO analysis, similar to the one that is customarily performed nowadays 
for deeply inelastic scattering, will be necessary. Moreover, by measuring the scale dependence at 
small x-Q, one can get a deeper insight into the structure of gluonic GPDs. 

There are several theoretical problems that should be answered in the future. 

The most crucial one is the magnitude of higher-twist corrections, especially, of the target 
mass corrections. Let us mention that the resummation of target mass effects from the twist- 


two operators performed in Ref. takes into account only a part of the entire tower of mass 
corrections, which are generated also by multi-particle operators in the off-forward Compton 
amplitude. They can be transformed into operators containing total derivatives by means of QCD 
equations of motion. Therefore, at any given nth order in the mass expansion, operators up to 
twist-n will contribute to the total correction. Therefore, the resummation cannot be handled in 
this situation, and it is sensible to address the problem order by order in the hadron mass. 

A deeper insight into the structure of perturbative series can be achieved by means of conformal 
symmetry 0,0. This formally requires to set the QCD /3-function to zero and to choose a special 
renormalization scheme, the so-called conformal scheme. With these simplihcations the Wilson 
coefficients of deeply inelastic scattering, known in next-to-next-to-leading order (NNLO), allow 
a numerical evaluation of the twist-two DVCS amplitude to NNLO in the valence quark region. 
Moreover, the corrections proportional to the /3-function have been calculated to the same order 
and can, therefore, be taken into account 0- 

If experiments can measure twist-three effects, one certainly needs their better theoretical 
understanding. Let us mention that the LO evolution of the antiquark-gluon-quark GPDs is 
governed by known two-particle kernels. Since the one-loop corrections to the short distance 
twist-two coefficient functions are important, one could anticipate that NLO effects can be as 
sizable in the twist-three amplitudes. This requires a detailed study of radiative corrections to 
the latter along the line similar to the computations done for the transversely polarized structure 
function g 2 of deeply inelastic scattering in Refs. [^, 0, 0, 0 . 


After these issues are addressed, their consequences and impact on the amplitudes have to be 
quantitatively analyzed in order to demonstrate the regions of the phase space, which are less 
sensitive to their presence and where an accurate determination of the parameters of the twist- 
two and -three GPDs can be performed. Of course, such studies are model-dependent. Thus, 
one should incorporate as much as possible information from perturbative QCD. This requires 
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modeling the GPDs at a low normalization point, e.g., Qo ~ 0.5 GeV, and their evolution upwards 
to experimental scales. 

In the most favorable situation, which means high enough required to discard power sup¬ 
pressed contributions and limit the analysis to the twist-two approximation in LO, the imaginary 
part of DVGS with real photon can give access to the GPDs on the diagonal x = ±,^ only, i.e., it 
becomes a function of a single scaling variable The only way one can measure it as 

a function of two variables, momentum fraction and skewedness, is in the deeply leptoproduction 
(virtual photoproduction) of a lepton pair, 

TVy* ^ N'-f* N'ii, 

since the generalized Bjorken variable and skewedness are independent in this case, ^ ^ f]- 

To conclude, experimental facilities having electron and positron beams is an ideal place to 
disentangle and study GPDs, and, thus, to extract fundamental information about the spin struc¬ 
ture of the nucleon. However, if one of them is not available, an accurate measurement with 
high momentum transfer probes will serve the purpose and will result into systematic tests of our 
understanding of the quark-gluon content of the nucleon. 
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